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SECTION  1 


SUMMARY 


1 . 1 BACKGROUND  AND  SCOPE 

In  1973  a concept  for  obtaining  higher  performance  from  two-stage  light  gas  launchers  was 
proposed  by  Physics  International  Company  (PI).  The  method  used  explosive  drivers  to  shock  heat 
a light-gas  working  medium  and  inject  it  ahead  of  a compressor  piston.  Calculations  indicated  that 
a significant  improvement  in  ballistic  performance  would  result,  making  possible  the  launching  of 
projectiles  three  to  four  times  more  massive  than  can  be  done  with  existing  launchers  of  the  same 
bore. 

From  April  1974  through  August  1975,  PI  conducted  a program  (Reference  1)  to  demonstrate  the 
concept  with  a launcher  of  2.54-cm  bore.  The  objective  of  that  program  was  to  launch  86-gm  projec- 
tiles to  velocities  in  the  range  of  3 to  6 km/sec.  However,  the  highest  velocity  obtained  in  eight 
tests  was  2.5  km/sec  and  the  compressor  chamber  and  barrel  of  the  launcher  were  damaged  by  erosion. 

The  Aerotherm  Division  of  Acurex  Corporation  investigated  the  PI  launcher  concept  for  the 
Defense  Nuclear  Agency.  Aerotherm's  objectives  were: 

• Evaluate  thermodynamic  and  gasdynamic  limitations  of  the  PI  hybrid  launcher 

• Recommend  configuration  and  materials  changes  that  promise  to  significantly  improve  per- 
formance 

This  report  presents  the  results  of  Aerotherm's  investigations. 

1.2  CONCLUSIONS 

The  observed  performance  shortfall  is  ascribed  to  the  particular  design  of  the  PI  apparatus, 
and  not  to  the  hybrid  concept. 

Specifically,  Aerotherm  has  concluded: 

• The  most  important  factor  in  performance  shortfall  observed  with  the  PI  launcher  was 
insufficient  working  medium  for  the  ballistic  cycle.  The  design  mass  of  hydrogen  was 
too  small,  and  a significant  quantity  failed  to  be  injected  into  the  compressor  tube. 

• The  maximum  attainable  kinetic  energy  of  the  75-mm  compressor  piston  is  marginal  for 
achieving  6 km/sec  muzzle  velocities 
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• The  junction  block  configurations  exceeded  the  maximum  angle(s)  for  favorable  injection 
of  the  driver  gas,  causing  reduction  of  the  driver  gas  injection  rate  for  large  driver 
to  compressor  diameter  ratios 

• Driver  injection  and  termination  dynamics  resulted  in  a significant  fraction  of  the 
driver  gases  being  trapped  between  the  termination  flanges  and  injection  ports 

• The  explosive  drivers  jetted,  which  may  have  led  to  contamination  of  the  driver  gas  and 
particle  erosion  of  internal  surfaces 

• Momentum  losses  through  the  boundary  layer  within  the  barrel  severely  limited  accelera- 
tion of  the  projectiles  beyond  300  barrel  diameters.  Pi's  design  calculations  excluded 
viscosity  effects,  and  because  of  this,  the  barrel  length  chosen  was  much  too  large  (576 
diameters) 

• Heat  losses  were  not  a major  cause  of  performance  shortfall;  they  could  account  for  no 
more  than  about  one-tenth  of  the  observed  reduction  of  muzzle  velocity 

• Computed  heat  loads  are  high  enough  to  cause  melting  of  the  steel  barrel's  inner  surface, 
particularly  near  the  chambrage  plane.  Melting  was  observed  in  the  experiments. 

• If  used  in  conjunction  with  properly  tailored  ballistic  cycles,  barrel  inserts/liners 
of  tungsten  (or  possibly  an  ablative  material  such  as  carbon)  can  protect  the  barrel 
from  destructive  melting 

t The  thermodynamic  description  of  the  hydrogen  driver  gas  must  account  for  molecular  co- 
volume and  vibration  as  well  as  dissociation 

! 
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SECTION  2 


INTRODUCTION 

The  PI  projectile  launcher  is  a hybrid  design  incorporating  a shock-heated  light-gas  working 
medium  with  a conventional  two-stage  ballistic  cycle.  In  principle,  the  preheated  working  medium 
leads  to  substantially  higher  reservoir  energy  densities  for  a given  volume  of  compressor  section. 
This  in  turn  leads  to  lower  pressure  gradients  throughout  the  launcher  for  any  ballistic  cycle,  in- 
cluding a constant  base  pressure  cycle.  As  a result,  substantially  higher  muzzle  velocities  for 
large  mass  to  diameter  projectiles  are  possible  within  the  pressure  limitations  of  launcher  design. 
With  the  PI  launcher,  however,  ballistic  performance,  gauged  by  projectile  muzzle  velocities,  was 
deficient. 

In  an  to  determine  the  causes  for  performance  shortfall,  the  following  possible  loss 

mechanism  'fied: 

jerly  operating  explosive  drivers 

• Incompatible  injection  geometry 

t Projecti le/barrel  friction 

• Momentum  losses  to  the  barrel  boundary  layer 

• Heat  losses  through  the  boundary  layer  and  consequent  erosion 

1 

: This  report  documents  the  Aerotherm  investigation  of  the  feasibility  of  the  PI  projectile 

launcher  as  follows: 

• Hydrogen  is  characterized  as  a Van  der  Waals  gas  subject  to  molecular  vibration  and 
dissociation 

t Baseline  predictions  define  ideal  launcher  performance 

• Parametric  results  show  performance  sensitivities  to  the  amount  of  working  medium  avail- 
able for  ballistic  cycles  and  to  various  other  important  ballistic  parameters 
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• The  loss  mechanisms  discussed  above  are  evaluated  for  impact  on  ballistic  performance 

• Material  response  is  determined  for  the  predicted  heat  loads 

On  the  basis  of  these  considerations,  design  modifications  to  achieve  desired  performance  levels  are 
recommended. 


SECTION  3 


r 


ESTIMATES  OF  PI  LAUNCHER  PERFORMANCE 

This  section  describes  a method  for  estimating  the  performance  of  the  PI  launcher  under 
idealized  conditions.  This  method  leads  to  reasonably  accurate  estimates  of  projectile  velocity  for 
representative  launcher  conditions  and  is  used  to  assess  the  sensitivity  of  muzzle  velocity  to  vari- 
ations of  key  ballistic  parameters.  These  idealized  performance  estimates  define  the  framework  within 
which  various  loss  mechanisms  can  be  evaluated. 

The  method  of  estimating  performance  is  based  on  a simple  thennodynamic  approach  to  describe 
the  shock  compression  and  injection  of  gas  by  the  explosive  drivers,  the  compression  of  the  gas  by 
the  first  stage  piston  to  a reservoir  state,  and  the  acceleration  of  the  projectile.  Throughout 
these  calculations,  an  accurate  equation  of  state  for  hydrogen  is  used  (see  Appendix  A). 

Performance  estimates  are  made  by  considering  the  operation  of  the  launcher  system.  The 
system  begins  at  an  ambient  state  characterized  by  a mass  of  gas  at  initial  pressure  and  temperature. 
The  explosive  drivers  are  then  actuated  to  process  the  driver  gas,  increasing  the  energy  of  the  gas 
by  an  amount  equal  to  the  work  done  by  the  collapsing  driver,  pAV.  The  driver  gas  is  considered  to 
be  completely  injected  into  the  compressor  section  and  the  compressor  piston  then  covers  the  driver 
gas  injection  ports.  The  compressor  piston,  having  been  accelerated  into  a vacuum  prior  to  reaching 
these  ports,  will  add  energy  to  the  gas  by  compression  as  it  decelerates  to  zero  velocity.  In  this 
model,  the  compressor  piston  decelerates  to  rest  before  the  projectile  is  allowed  to  accelerate.  The 
final  compressed  state  is  then  used  to  estimate  projectile  velocity  with  a preburned-propel 1 ant  bal- 
listic model . 

In  reality,  the  projectile  begins  to  accelerate  before  the  compressor  piston  decelerates  to 
rest.  Consequently,  the  method  used  here  somewhat  overestimates  the  reservoir  pressures,  temperatures, 
and  sound  speeds.  However,  for  a sufficiently  long  launch  tube,  the  estimated  muzzle  velocity  using 
this  method  is  accurate.* 

The  various  steps  in  the  performance  calculation  are  now  described  in  detail. 


★ 

Future  performance  estimates  should  be  made  using  a 1-1/2-D  finite  difference  code  incorporating 
the  more  accurate  equation  of  state  for  hydrogen  used  in  this  report. 


3.1 


SHOCKED  GAS  STATE  IN  THE  EXPLOSIVE  DRIVERS 


The  explosive  drivers  can  be  modeled  as  a constant  velocity  piston  (at  the  detonation  velocity 
of  the  explosive)  which  drives  a shock  wave  into  the  driver  gas.  Knowing  the  initial  driver  gas 
pressure  and  temperature  and  by  measuring  the  shock  velocity,  the  gas  conditions  behind  the  shock 
can  be  calculated  by  solving  the  conservation  equations  and  the  equation  of  state  of  the  gas. 


The  conservation  equations  across  the  shock,  expressed  in  shock  coordinates  are; 


PlVi  P2V2 


Po  + 


P2  V2^  + p^  v.| 


P2  '1 2 


P,  V 2 


(1) 


where  p = pressure 
p = density 

V = particle  velocity  in  shock  coordinates 

6 = ^ p/p'~*^’  P^i^aneter  introduced  for  convenience  to  preserve  the  ideal  form  of  the 

shock  jump  relations  for  a nonideal  gas) 
and  e = specific  internal  energy  (see  Appendix  A) 


The  equation  of  state  is: 


p = p(T,  i,r.a ,b,  •^)  (see  Appendix  A) 

where  T = temperature 

a = degree  of  dissociation 
a = Van  der  Waal's  attraction  constant 
b = covolume  constant 

9^  = characteristic  vibrational  temperature 


These  equations  are  solved  by  an  iterative  process  to  obtain  P2 , P2,  62,  and  v^,  and  trans- 
formed back  to  laboratory  coordinates  by 


u 


2 


s 


where  U2  = particle  velocity  behind  the  shock  in  laboratory  coordinates 
s = observed  shock  velocity  in  laboratory  coordinates 
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The  shocked  gas  states  for  various  shock  wave  velocities  (6.0  km/sec  ■.  s 10.0  km/sec)  into 
hydrogen  initially  at  300“K  and  350  psia  are  sumniarized  in  Table  1.  Intermolecular  attraction  and 
dissociation  are  negligible,  but  covolume  and  molecular  vibration  play  an  important  role.  Pressure 


behind  the  shock  is  relatively  insensitive  to  covolume  forces,  but  flow  Mach  number  and  density  ratio 
across  the  shock  are  significantly  affected  by  covolume  forces.  The  calculated  variation  of  flow 
Mach  number  and  density  ratio  with  shock  velocity  are  shown  in  Figure  1 with  and  without  covolume 
forces.  Calculated  particle  velocity  and  shocked  pressure  variation  with  shock  velocity  are  shown 
in  Figure  2. 

When  the  explosive  driver  is  operated  within  its  range  of  ideal  performance,  the  particle 
velocity  behind  the  shock  is  equal  to  the  detonation  velocity  of  the  explosive  and  conditions  behind 
the  shock  are  constant.  PI  used  nitromethane  explosive  (D  = 6.25  km/sec).  For  hydrogen  initially 
at  300°K  and  350  psia,  the  shock  velocity  should  be  7.8  km/seci  thus,  in  subsequent  launcher  perfor- 
mance estimates,  the  explosive  drivers  are  assumed  to  operate  ideally  and  shocked  gas  conditions  cor- 
responding to  a driver's  shock  velocity  of  7.8  km/sec  are  used. 

3.2  CALCULATION  OF  RESERVOIR  CONDITIONS 

The  initial  state  of  the  PI  launcher  system  is  taken  as 


where  p^  = initial  driver  pressure 

Pj  = initial  compressor  pressure 
T.|  = ambient  temperature 
Up  = projectile  velocity 

u^^j  = compressor  piston  velocity  at  injection  port 

From  the  initial  temperature  and  loading  pressures  the  mass  of  gas,  G,  can  be  determined.  The 
initial  energy  of  the  system  is  then  given  by: 

E.  = 5/2  RT^  X G (2) 


where  R = gas  constant  of  molecular  hydrogen 


TABLE  1.  CONDITIONS  BEHIND  A SHOCK  IN  HYDROGEN  INITIALLY  AT  350  PSI  AND  300°K 


Full 

Equati on 

of  State 

s 

(km/sec) 

T2 

(°K) 

P2 

(gm/cm’) 

P2 

kb 

“2 

'^adiab 

^isent 

"2 

(km/ sec) 

<^2 

(km/sec) 

? 

Pj/P, 

9.0 

2699 

0.0104 

1.31 

0.0009 

1.410 

1.497 

4.35 

7.31  1. 

58 

5.31 

8.5 

2467 

0.0101 

1.16 

0.0004 

1.415 

1.497 

4.15 

6.86  l.f 

55 

5.20 

8.0 

2239 

0.0099 

1.03 

0.0001 

1.420 

1.501 

3.96 

6.42  !.( 

52 

7.5 

2019 

0.0096 

0.90 

<0.0001 

1.424 

1.503 

3.75 

5.98  1. 

59 

4.95 

7.0 

1810 

0.0094 

0.78 

1.427 

1.509 

3.55 

5.55  1. 

56 

4.81 

6.5 

1610 

0.0091 

0.67 

1.431 

1.514 

3.35 

5.11  1. 

53 

4.67 

6.0 

1422 

0.0088 

0.57 

1.433 

1.519 

3.15 

4.67  1.^ 

18 

4.51 

Covolume  Suppressed 

9.0 

2804 

0.0116 

1.34 

0.0012 

1.358 

1.316 

3.90 

7.49  l.< 

32 

5.96 

8.5 

2561 

0.0113 

1.19 

0.0005 

1.364 

1.321 

3.74 

7.04  1. 

38 

5.81 

8.0 

2324 

0.0110 

1.05 

0.0002 

1 . 369 

1.327 

3.57 

6.59  1.! 

35 

5.66 

7.5 

2096 

0.0107 

0.92 

<0.0001 

1.375 

1.334 

3.40 

6.14  1. 

31 

5.50 

7.0 

1878 

0.0104 

0.80 

1.380 

1.343 

3.22 

5.69  1. 

77 

5.31 

6.5 

1671 

0.0100 

0.69 

1.385 

1.352 

3.05 

5.24  1. 

72 

5.12 

6.0 

1476 

0.0097 

0.58 

1.388 

1.362 

2.88 

4.79  l.( 

56 

4.94 

12 


When  the  explosive  drivers  are  actuated,  the  volume  of  the  system  changes  and  the  explosively 


formed  "pistons"  do  work  at  constant  pressure.  After  driver  operation  is  complete,  the  system  has 
a total  energy  of: 

Eg  = E.  + PgAV  (3) 

where  AV  = change  in  volume  for  both  drivers  resulting  from  pressure  tube  collapse 
Pg  = pressure  behind  the  driver  shock 

All  gas  processed  by  the  explosive  drivers  is  now  assumed  to  be  carried  by  its  own  momentum  into  the 
compressor  section  and  trapped  there  when  the  compressor  piston  just  covers  the  injection  ports. 

The  compressor  piston,  separated  from  the  rest  of  the  system  by  a diaphragm,  is  accelerated 
into  a vacuum.  In  moving  from  the  diaphragm  position  just  upstream  of  the  injection  ports  to  the 
point  where  it  just  covers  the  injection  ports,  trapping  the  injected  gas,  the  compressor  piston 
does  negligible  work  on  the  system.  System  energy  is  still  accurately  given  by  Equation  (3). 

At  this  point,  the  system  is  considered  to  be  in  a state  shown  by: 


2 

where  e,  = -?- 

o u 

p = density  of  the  gas,  G/V.  • 

J .«»  * '*J 

' volume  of  the  compressor  section  at  injection  port  closure) 

These  two  state  variables  are  sufficient  to  determine  all  other  state  variables  in  the  system  using 
Figures  A-1  thrcugh  A-5  in  Appendix  A,  or  numerically  using  relations  (A-1)  to  (A-7),  also  in  Appendix  A. 

Compression  is  assumed  to  take  place  isentropical ly . and  work  is  done  on  the  system  by  the 
compressor  piston  as  it  decelerates  and  gives  up  its  kinetic  energy  to  the  gas.  Compression  begins 
at  piston  velocities  of  3,500  fps  (1.1  km/sec)  or  less,  so  that  the  kinetic  energy  of  the  gas  is 
always  much  less  than  its  internal  energy.  Compression  can  be  modeled  as  a series  of  equilibrium 
states . 
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Pressures  and  densities  at  constant  entropy  are  determined  graphically  from  Figures  A-1  to  A-5, 


or  numerically  from  relations  (A-1)  through  (A-7). 

As  the  compressor  piston  decelerates,  the  gas  density  increases  as  the  volume  of  the  compressor 
section  decreases.  The  relation  between  gas  density  and  compressor  section  geometry  is 


where  A = cross  sectional  area  of  compressor  section 

X = distance  measured  ^rom  downsti^eam  edge  of  injection  port 

The  pressures  and  densities  at  constant  entropy  are  fitted  to  a power  relation  of  the  form* 


(4) 


p = np^ 


(5) 


Using  (4)  and  (5),  the  deceleration  of  the  piston  is  calculated  by  integrating  the  expression 


pdx  = (1/2) 


pist 


where  p^  = pressure  when  compressor  piston  stops 

Pj  = pressure  after  injection  but  before  compression 
Mpist  = compressor  piston  mass 

Performing  the  indicated  integration  to  the  position  where  the  piston  comes  to  a stop: 


(6) 


V.  . 
inj 


(y-1)  M, 


pist 


+ V, 


(1-y) 


n J_ 

1-Y 


2n^ 


inj 


PS 


where  ps  indicates  "piston  stop." 


(7) 


This  equation  gives  the  distance  from  the  downstream  edge  of  the  injection  port  to  the  point 
at  which  the  piston  stops,  assuming  full  injection,  no  losses,  and  no  projectile  motion.  However, 
in  reality,  complete  injection  never  happens,  losses  occur,  and  the  projectile  begins  to  move  before 
the  compressor  piston  stops.  As  a result,  the  piston  will  stop  further  downstream  than  indicated  by 
(7). 

★ 

This  equation  gives  an  accurate  representation  of  pressure  variation  with  density  for  calculating 
the  deceleration  of  the  compressor  piston.  This  is  the  only  place  it  is  used  in  this  procedure. 
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Once  Xpj  is  determined,  reservoir  density,  is  determined  from  (4).  Since  compression  was 
assumed  to  be  isentropic,  (where  S = entropy)  and  the  reservoir  state  is  completely  determined 

from  the  two  state  variables  p^,  S^,  and  Equations  (A-1)  through  (A-7).  Alternately,  the  reservoir  state 
can  be  determined  graphically  from  Figures  (A-1)  through  (A-5). 

Dimensions  and  volumes  for  various  configurations  used  by  PI  are  given  in  Figure  3 and  Table  2. 
Calculated  reservoir  states  are  given  in  Table  3 for  several  configurations  and  initial  conditions 
of  the  PI  launcher. 

Since  the  hydrogen  is  partially  dissociated  and  influenced  by  covolume  forces  and  molecular 
vibration,  the  sound  speeds  and  isentropic  indices  were  computed  numerically  from  relations  (A-1), 

(A-2) , and  (A-5)  using 


AS  = 0 
AT  = 1°K 


( s.\  r^i 

\ p / 


AS  = 0 
AT  = 1°K 


where  C = sound  speed 

Yisent  ” isentropic  index 
AS  = change  in  entropy 

Reservoir  temperature,  T^,  is  shown  in  Figure  4 as  a function  of  reservoir  density,  for 
■Tcveral  pressures  and  scaled  sound  speeds,  C^/y^.  Figure  4 is  used  later  in  estimating  heat  loss 
from  the  reservoir. 

3.3  CALCULATION  OF  BALLISTIC  CYCLE  PERFORMANCE 


The  calculation  of  projectile  velocity  is  made  by  isentropically  expanding  the  gas  from  the 
reservoir  state  using  the  preburned-propellant  gun  model  and  curves  of  Seigel  (Reference  2).  In 
this  model,  the  compressor  piston  is  assumed  to  have  decelerated  to  rest  permanently,  prior  to  pro- 
jectile acceleration,  hence  the  reservoir  geometry  is  constant.  In  using  Seigel 's  preburned-propel- 
lant gun  results,  the  reservoir  is  characterized  by  pressure,  p^ ; sound  speed,  C^;  gas  mass,  G; 
isentropic  index,  y^;  and  compressor  to  barrel  diameter  ratio,  02^.  The  launch  tube  has  length  L and 
the  projectile  is  characterized  by  its  mass,  Mp. 

Typical  performance  curves  are  shown  in  Figures  5 to  7.  These  curves  were  derived  by  computer 

solution  of  the  characteristic  equations  for  various  gas  to  projectile  mass  ratios,  G/Mp,  and  the 

reservoir  to  launch  tube  diameter  ratio,  D,,.  The  curves  are  displayed  with  the  dimensionless  pro- 

_ P4AL 

jectile  velocity  u = plotted  against  the  ballistic  parameter  x = ' 

" p'“4 
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Figure  3.  PI  hybrid  two  stage  launcher  geometry  (see  Table  2). 
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TABLE  2.  DIMENSIONS  AND  VOLUMES  FOR  PI  LAUNCHER  TESTS  (SEE  FIGURE  3) 


Test  2 

Test  3 

Test  4 

Test  5 

Test  6 

Test  7 

Test  8 

A 

417.83 

417.83 

417.83 

417.83 

1 B 

1 

20.32 

20.32 

20.32 

38.74 

C 394.97 

394.97 

394.97 

394.97 

394.97 

394.97 

394.97 

394.97 

D 35.56 

35.56 

35.56 

35.56 

35.56 

35.56 

35.56 

45.72 

E 

240.94 

240.94 

240.94 

240.94 

F 18.03 

18.03 

18.03 

18.03 

18.03 

18.03 

18.03 

18.03 

G 

19.69 

19.69 

19.69 

H 

5.08 

5.08 

5.08 

J 

11.43 

11.43 

11.43 

L 22.86 

22.86 

50.80 

50.80 

50.80 

50.80 

57.15 

57.15 

V 

400.05 

400.05 

400.05 

0-  (degrees) 

45° 

45° 

45° 

45° 

O 

uo 

45° 

24° 

Volume 

(cm’) 

Driver  pressure 
tubes  to  flange 

25,963 

25,963 

25,963 

12,982 

Driver  flange  to 
compressor  port 

1 ,287 

1 ,287 

1 ,287 

1,227 

Compressor 
between  diaphragm 
and  projectile 

13,304 

13,304 

13,304 

13,780 

Compressor  after 
port  closure 

12,352 

12,352 

12,352 

11  ,710 

AV  due  to 
driver  operation 

25,242 

25,242 

25,242 

12,621 

9 


Reservoir  temperature  (°K) 


- lUiSC 


These  curves  are  available  for  = 1-1,  1-4,  and  1.67.  In  the  range  of  conditions  in 

which  the  PI  launcher  operates,  the  reservoir  gas  is  vibrational ly  excited,  partially  dissociated, 
and  is  in  a density  range  where  covolume  forces  are  significant.  These  reservoir  conditions  must  be 
converted  into  equivalent  conditions  compatible  with  Siegel's  curves. 

Glass  (Reference  3)  has  shown  that  muzzle  velocities  are  nearly  invariant  for  reservoir  gases 
of  the  same  scaled  sound  speed,  C^/y^  where  the  condition  Up  < 2 is  satisfied.  For  conditions 

typical  of  the  PI  launcher,  this  condition  is  always  met.  : 

j 

Thus,  for  a reservoir  with  a calculated  and  y^,  an  equivalent  reservoir  with  and  y^  = I 

1.1,  1.4,  or  1.67  can  be  specified  where  C^/y^  = C^/y'^.  i 

The  Seigel  curves  assume  y^  is  constant  throughout  the  projectile  acceleration.  When  covolume 
forces  are  not  present,  this  assumption  is  valid  over  the  range  of  operating  conditions  characteris- 
tic of  the  PI  launcher  (Reference  4).  When  covolume  forces  are  important,  this  assumption  leads  to 
slightly  lower  calculated  muzzle  velocities. 

Based  on  the  calculation  procedure  illustrated  in  Appendix  B,  performance  estimates  for  the 
initial  conditions  of  Table  3 are  presented  in  Table  4.  Some  conclusions  are  worth  noting: 

t Evacuating  the  compressor  section  results  in  a smaller  G/Mp  but  a higher  reservoir 
specific  energy.  Projectile  velocity  is  only  slightly  higher. 

• Substantial  performance  increases  result  if  piston  energy  is  increased 

• Optimum  performance  is  achieved  if  the  compressor  section  is  evacuated,  piston  energy 
is  maximized,  and  if  the  G/Mp  ratio  is  greater  than  1 

The  upper  limit  on  muzzle  velocity  for  the  present  configuration  of  the  PI  launcher  can  be  es- 
timated by  this  method.  As  jioted  previously,  this  estimate  assumes  ideal  driver  operation,  full  in- 
jection, and  no  loss  mechanisms.  To  estimate  the  upper  limit,  the  following  initial  conditions  are 
used: 

• Initial  driver  pressure  is  doubled  to  700  psia  to  increase  the  G/Mp  ratio 

• The  compressor  chamber  is  evacuated  to  maximize  reservoir  specific  energy 

t The  75  mm  gun  is  operated  at  its  performance  limit  to  maximize  compressor  piston  energy 

The  state  variables  and  conditions  which  are  achieved  through  these  modifications  are  pre- 
sented in  Tables  3 and  4 (Table  Index  T8).  The  calculated  muzzle  velocity  for  the  86  gram  projectile 


TABLE  3.  CALCULATED”  DRIVER,  COMPRESSOR,  AND  RESERVOIR  STATES  FOR  SEVERAL  DRIVER/COMPRESSOR  CONDITIONS 


ni tial 
'•essure 
(Psi) 

Shock 

Pressure 

(kilobars) 

Piston 

Mass 

(cms) 

Piston 

Velocity 

(ktn/sec) 

Initial 

Pressure 

(psi) 

Max  piston  stop 
from  transition 
(cm) 

*’4 

(kilcbars) 

^4 

(“K) 

^4 

(gm/cm* 

350  0.974  2724  0.91  350  101.35  2.98  3715  0.0159 

Test  7 

350 

0.953 

2724 

0.91 

350 

84.58 

2.89 

3715 

0.0187 

Test  7, 

covolume 

suppressed 

350 

0.974 

2724 

0.91 

3.7 

81.79 

2.94 

4785 

0.0124 

350 

0.974 

3450 

0.91 

3.7 

63.25 

4.07 

5050 

0.0156 

350 

0.974 

2724 

0.91 

350 

52.32 

3.88 

3613 

0.0201 

Test  8 

1 driver 

700 

1.99 

4540 

0.91 

350 

118.36 

5.20 

4155 

0.0226 

700 

1.99 

4540 

0.91 

3.7 

105.92 

5.05 

4803 

0.0196 

700 

1.99 

4540 

1.07 

3.7 

81.03 

7.40 

5120 

0.0250 

700 

1.99 

4540 

0.76 

^.7 

134.87 

3.61 

4534 

0.0156 

700 

1.99  . 

4540 

1.22 

J.7 

64.01 

10.47 

5425 

0.0308 

TABLE  4.  CALCULATED  LAUNCHER  PERFORMANCE  FOR  SEVERAL  DPIVER/COMPRESSOR 
CONDITIONS  (RESERVOIR  CONDITIONS  OF  TABLE  3) 


Driver 

Ini tial 
Pressure 
(psi) 

Piston 

Mass 

(gms) 

Compressor 

Piston  Initial 
Velocity  Pressure 
(km/sec)  (psi) 

Reservoi r 
P4 

(Kilobars)  (km/sec) 

Launch  Cycle 

*^P  “^P 

(gms)  (km/sec) 

Table 

Index 

350 

2724 

0.91 

350 

2.98 

3.51 

86 

0.96 

3.70 

T1 

Like  Test 

7*^ 

350 

2724 

0.91 

350 

2.89 

3.52 

86 

0.96 

3.67 

T2 

Covolume 

suppressed 

350 

2724 

0.91 

3.7 

2.94 

4.11 

86 

0.62 

3.80 

T3 

350 

3450 

0.91 

3.7 

4.07 

4.28 

86 

0.62 

4.28 

T4 

350 

2724 

0.91 

350 

3.88 

3.45 

87.4 

0.64 

3.74 

T5 

Like  Test 

8^ 

700 

4540 

0.91 

350 

5.20 

3.72 

86 

1.57 

4.65 

T6 

700 

4540 

0.91 

3.7 

5.05 

4.14 

86 

1.23 

4.81 

T7 

700 

4540 

1.07 

3.7 

7.40 

4.32 

86 

1.23 

5.39 

T8 

700 

4540 

0.76 

3.7 

3.61 

4.00 

86 

1.23 

4.28 

T9 

700 

4540 

1 .22 

3.7 

10.47 

4.46 

86 

1.23 

6.12 

TIO 

Calculation  procedure  is  described  in  Appendix  B 
For  L/D  = 300 

Experimental  ve’locities  were  about  70  percent  of  these  ideal  results 


s: 


is  5.39  km/sec.  Thus  in  the  present  configuration  (75  mm  gun,  2 drivers),  the  upper  performance 
limit  under  ideal  conditions  appears  to  be  about  5.5  km/sec  (18,000  fps). 

Reconfiguration  of  the  system,  such  as  addition  of  another  driver  and  use  of  a larger  first 
stage,  is  necessary  to  increase  the  performance  limit  to  the  design  goal  of  6.1  km/sec.  In  fact,  the 
performance  limit  calculated  under  the  above  idealized  conditions  should  be  somewhat  greater  than 
6.1  km/sec  to  provide  some  margin  for  loss  mechanisms  such  as  incomplete  gas  injection,  reservoir 
heat  loss,  and  boundary  layer  effects. 

3.4  SENSITIVITY  OF  MUZZLE  VELOCITY  TO  RESERVOIR  AND  LAUNCHER  PARAMETERS 

The  sensitivity  of  muzzle  velocity  to  the  following  parameters  has  been  evaluated: 

t Barrel  length  to  diameter  ratio,  L/D 

• Hydroaen  gas  to  projectile  mass  ratio,  G/Mp 

• Compressor  piston  energy 

• Reservoir  pressure,  p^  and  scaled  sound  speed 

Figures  8 and  9 show  the  variation  of  muzzle  velocity  with  barrel  L/O  for  conditions  typical 
of  the  PI  tests.  As  will  be  discussed  in  Section  4,  projectile  velocity  is  not  expected  to  increase 
significantly  beyond  a barrel  L/D  of  about  300  because  of  boundary  layer  effects.  The  performance 
predicted  by  PI,  however,  depends  on  the  velocity  calculated  for  a barrel  L/D  of  576.  For  conditions 
representative  of  the  PI  tests,  the  last  half  of  the  barrel  would  contribute  an  additional  15  to  20 
percent  to  projectile  velocity,  if  boundary  layer  effects  did  not  exist.  It  is  therefore  important 
in  future  work  to  use  barrel  lengths  no  greater  than  about  300  L/D's  in  preliminary  performance 
estimates.*  In  assessing  the  sensitivity  of  muzzle  velocity  to  varying  G/M^,  p^,  and  C^/y^,  the 
barrel  length  is  henceforth  considered  fixed  at  300  L/D's. 

Figures  10  and  11  cover  the  range  of  p^  and  C^/y^  of  the  PI  tests;  they  show  that  muzzle  velocity 
is  quite  sensitive  to  the  G/Mp  ratio  for  the  range  of  design  G/Mp  values  in  the  PI  tests  (0.62 
G/Mp  < 0.92  without  any  losses).  From  these  results  it  can  be  concluded  that  any  loss  of  driver  gas 
during  driver  operation,  injection,  compression,  or  initial  projectile  acceleration  in  the  PI  tests 
would  seriously  degrade  anticipated  performance.  A G/Mp  ratio  of  1.5  to  2.0  would  seem  to  be  a good 
design  goal  since  small  gas  losses  could  be  incurred  without  seriously  degrading  performance. 


However,  no  performance  penalty  should  result  from  experimenting  with  a large  L/D  launcher  to  derive 
experimental  data  on  the  benefits  of  incremental  lengths.  Similarly,  larger  L/D's  could  be  con- 
sidered when  boundary  layer  effects  are  adequately  accounted  for  in  the  evaluation  procedure. 
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Figure  9.  Performance  sensitivity  to  barrel  dimensions 
for  various  ratios  of  gas  mass  to  projectile 
mass,  p.  - 10  kb. 
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Performance  sensitivity  to  ratio  of  gas  mass 
to  projectile  mass  for  representative  scaled 


reservoir  sound  speeds,  P/,  = 5 kb. 
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Figure  11.  Performance  sensitivity  to  ratio  of  gas  mass 
to  projectile  mass  for  representative  scaled 
reservoir  sound  speeds,  p.  = 10  kb. 


Further  increases  in  the  G/Mp  ratio  would  not  lead  to  substantially  better  performance,  but  would 
provide  a larger  margin  against  the  inevitable  gas  losses. 

Another  important  design  goal  would  be  to  maximize  the  compressor  piston  energy,  since  piston 
energy  has  a strong  effect  on  muzzle  velocity.  (For  example,  piston  energy  controls  final  reservoir 
density  and  adds  substantial  energy  to  the  reservoir.)  Figure  12  shows  that  under  typical  operating 
conditions,  a 40-percent  increase  in  piston  velocity  results  in  a 26-percent  increase  in  muzzle 
velocity. 

The  method  used  herein  is  accurate  for  estimating  muzzle  velocities  and  assessing  the  sensi- 
tivity of  performance  to  variations  in  G/Mp  ratio,  scaled  sound  speed,  C^/y^,  etc.  In  these  calcula- 
tions, it  is  assumed  that  the  projectile  is  restrained  until  the  compressor  piston  has  come  to  rest. 
Actually,  the  projectile  is  well  into  the  launch  cycle  before  the  compressor  piston  comes  to  rest, 
so  that  peak  reservoir  pressures  and  sound  speeds  are  considerably  lower  than  those  calculated.  The 
calculated  pressure  levels  shown  in  Figure  13  are  thus  unrealistically  high,  but  the  figure  does 
illustrate  the  sensitivity  of  muzzle  velocity  to  peak  pressure.  If  the  G/Mp  ratio  is  adequate  (i.e., 
G/Mp  >1),  the  failure  to  develop  full  design  pressures  will  not  seriously  affect  muzzle  velocity, 
provided  that  the  scaled  sound  speed  remains  high.  At  G/Mp  levels  much  less  than  1.0,  as  in  some  of 
the  PI  experiments,  performance  is  quite  sensitive  to  reservoir  pressure  levels. 

Figures  14  and  15  illustrate  the  sensitivity  of  muzzle  velocity  to  scaled  sound  speed,  C^/y^ 
in  the  reservoir.  The  combination  of  high  sound  speed,  and  low  isentropic  index,  y^,  leads  to 
the  highest  muzzle  velocities.  The  scaled  reservoir  sound  speed  is  an  important  parameter.  It  can 
be  substantially  lowered  by  heat  transfer  losses  during  the  launch  cycle  and  lowered  further  if  heat 
transfer  is  severe  enough  to  trigger  erosion  and  ablation. 

In  a gas  where  covolume  effects  are  important,  both  sound  speed,  C^,  and  isentropic  index, 
i increase  rapidly  with  further  increases  in  density.  However,  the  scaled  sound  speed,  C^/>^  in- 
creases gradually  and  closely  follows  the  rise  in  internal  energy,  e^. 

In  summary,  these  studies  show  that  projectile  velocity  increases  with  barrel  length,  gas-to- 
mass  ratio,  compressor  piston  velocity,  reservoir  pressure,  and  scaled  speed  of  sound  in  the  reser- 
voir. The  L/D  effect  is  small  beyond  300  even  when  boundary  layer  effects  are  ignored.  On  the  other 
hand,  performance  is  quite  sensitive  to  gas-to-mass  ratio  below  unity,  which  was  typical  of  the  PI 
experiments . 
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Figure  14.  Performance  sensitivity  to  scaled  reservoir  sound 
speed  for  various  ratios  of  gas  mass  to  projectile 
mass,  p^  = 5 kb. 
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Figure  15.  Performance  sensitivity  to  scaled  reservoir  sound 
speed  for  various  ratios  of  gas  mass  to  projectile 
mass,  p,  = 10  kb. 


Substantial  gains  in  performance  can  be  expected  by  utilizing  the  maximum  capability  of  the 
75  mm  gun  to  benefit  from  increased  peak  chamber  pressure  and  scaled  speed  of  sound.  However,  the 


F 


capability  of  the  75  mm  gun  is  marginal  in  terms  of  achieving 
for  an  86  gm  projectile  in  a 2.54-cm  (1-inch)  diameter  launch 


the  desired  6.1  km/sec  muzzle  velocity 
tube. 
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SECTION  4 


FACTORS  DEGRADING  PERFORMANCE 

This  section  will  deal  with  the  difficulties  in  obtaining  the  desired  chamber  conditions  in 
the  launcher,  based  on  the  driver  and  injection  processes.  Losses  accrued  during  the  ballistic  cycle 
after  completion  of  injection  will  also  be  evaluated. 

4.1  GAS  CONDITIONS  GENERATED  BY  THE  EXPLOSIVE  DRIVERS 

In  the  ideal  mode  of  explosive  driver  operation  shown  in  Figure  16  (no  jetting,  no  gas  loss, 
no  radial  expansion),  the  particle  velocity  behind  the  shock  is  equal  to  the  explosive  driver  detona- 
tion velocity.  The  shockec  as  density  and  pressure  remain  constant  behind  the  shock. 

In  practice,  this  mode  of  operation  can  be  maintained  up  to  approximately  10  pressure  tube 
diameters  of  shocked  gas.  Thereafter,  the  momentum  of  boundary  layer  gases  is  sufficient  to  prevent 
complete  collapse  of  the  pressure  tube.  The  system  then  begins  to  lose  driver  gas,  with  an  accompany- 
ing decrease  in  shock  velocity  after  a suitable  communication  time. 

During  ideal  operation,  the  volume  of  shocked  gas  compressed  by  the  shock  wave  is  equal  to 
the  total  volume  processed  by  the  shock,  divided  by  the  compression  ratio  across  the  shock.  Also, 
to  properly  determine  the  position  of  the  shock  relative  to  the  detonation  wave,  the  volume  of  the 
collapse  cone  must  be  included.  This  is  true  because  with  a typical  collapse  half-angle  of  9°,  the 
volume  of  the  collapse  cone  is  about  10  percent  of  the  total  volume  for  the  maximum  length  of  shocked 
gas  in  ideal  drivers. 

It  can  be  shown  that  the  distance  between  the  shock  and  the  detonation  wave,  should  be 


Ax 


p^i  6 tan  o 


where 


x^  = distance  shock  has  travelled  in  unshocked  gas 
p^i  = compression  ratio 


(8) 
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Dj.  = pressure  tube  inside  diameter 
0 = collapse  half-angle 

A driver  will  jet  when  the  collapse  process  is  overdriven.  Under  these  conditions,  the  maxi- 
mum amount  of  jet  material  possible  is  determined  from  the  classical  jetting  formula  (Reference  5): 


"li  = 1 - cos  o 
'"a  ' ^ 


where  m^  = mass  of  jet 

= mass  of  liner  or  pressure  tube 

In  practice,  however,  the  maximum  jet  mass  is  never  realized  since  driver  gas  pressure,  ir- 
reversible plastic  working  of  the  liner  during  collapse,  and  liner  compressibility  effects  all  tend 
to  suppress  jetting. 

For  a nitromethane  driver  (detonation  velocity  D = 6.25  km/sec)  operating  in  its  ideal  range, 
the  calculated  shock  velocity  is  7.8  km/sec.  The  calculated  compression  ratio  across  the  shock  is 
5.03  for  hydrogen  initially  at  300°K  and  350  psia.  Examination  of  the  PI  data  for  the  drivers  used 
in  Tests  5 to  8 (Figures  17  to  20)  shows  abnormally  high  shock  velocities  ranging  from  8.6  to  9.0 
ktn/sec,  characteristic  of  severely  jetting  drivers. 

In  Test  7 for  example,  the  observed  shock  velocity  is  8.7  km/sec  and  the  distance  between  the 
shock  and  detonation  wave  (Figure  19)  leads  to  an  apparent  compression  ratio  of  3.48  (Equation  (8) 
for  a = 9°),  if  shocked  gas  conditions  are  assumed  uniform.  For  a shock  velocity  of  8.7  km/sec, 
however,  a compression  ratio  of  5.2  is  calculated  (Figure  1).  Therefore,  shocked  gas  conditions 
for  this  driver  must  be  very  nonuniform. 

The  observed  shock  velocities  are  constant  over  the  length  of  the  drivers  and  are  reproduc- 
ible from  driver  to  driver.  The  most  likely  explanation  of  this  type  of  reproducible  nonideal  per- 
formance is  that  the  collapsing  pressure  tube  is  overdriven  and  jets  a spray  of  metal  particles 
into  the  driver  gas.  At  the  pressures  and  temperatures  involved,  the  jet  particles  would  be  solid 
or  liquid  and  their  volume  would  be  negligible.  However,  the  particles  are  jetted  into  the  driver 
gas  at  about  twice  detonation  velocity  (i.e.,  about  12.5  km/sec)  in  laboratory  coordinates.  As 
these  particles  are  decelerated  by  the  shocked  gas,  they  transfer  momentum  to  the  gas  and  thereby 
can  support  pressure  and  density  gradients. 

If  it  is  assumed  that  this  expanding  cloud  of  steel  particles  "snowplows"  and  forms  an  im- 
permeable piston,  then  this  particle  piston  would  have  to  support  a pressure  differential  of  1.2  kb 
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Figure  19.  Explosive  driver  performance  for  PI,  Test  Number  7. 
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to  produce  the  observed  shock  velocity  of  8.7  kni/sec.  By  conservation  of  momentum,  the  flux  of  jet 
particles  would  be  80  percent  of  the  theoretical  maximum  established  by  Equation  (9)  for  o = 9° 
and  mj^  = 16.4  gms/us. 

However,  it  is  more  likely  that  the  jet  particles  would  be  distributed  throughout  the  shocked 
gas.  Initially  jet  particles  would  penetrate  the  gas  ahead  of  the  shock,  then  as  the  shock  moves 
downstream,  it  would  engulf  the  original  particles.  Meanwhile,  the  constant  flux  of  new  particles 
would  drag  decelerate  to  local  flow  velocity  and  travel  with  the  gas  behind  the  shock.  This  process 
is  illustrated  schematically  in  Figure  21.  While  the  particles  can  be  viewed  as  decelerating  in  the 
gas  by  drag,  the  gas  can  be  viewed  as  flowing  through  a constant  area  duct  with  friction  (the  sum 
of  particle  cross  sections  is  negligible  compared  to  the  pressure  tube  cross  section,  and  the  total 
particle  drag  can  be  considered  as  an  average  friction  factor).  Thus,  pressure  and  density  would 
decrease  from  a maximum  behind  the  shock  and  the  flow  Mach  number  relative  to  the  shock  would  in- 
crease towards  unity.  The  pressure  drop  through  the  gas  would  reflect  the  jet  momentum  input.  The 
density  drop  through  the  gas  would  explain  the  apparent  increase  in  length  of  the  shocked  gas. 

Assuming  a particle  diameter  of  1.0  mm,  the  maximum  possible  flux  of  jet  material  would  pro- 
duce 24  particles  per  microsecond  (ps).  This  production  rate  seems  too  low  to  account  for  the  ob- 
served uniformity  and  reproducibility  of  the  shock  trajectories  in  the  PI  tests.  An  assumed  particle 
diameter  of  0.1  mn  would  result  in  a production  of  24,000  particles  per  ps  (which  seems  to  be  a more 
reasonable  value)  for  the  theoretical  maximum  flux  of  jet  material. 

If  these  solid  or  liquid  particles  were  roughly  spherical,  the  drag  coefficient  would  be 
about  0.4  and  constant  (10’  < Reynolds  number  < 10’).  A typical  deceleration  time  from  12.5  km/sec 
to  7 km/sec  in  the  shocked  hydrogen  is  250  ps.  A deceleration  time  of  250  ps  is  consistent  with 
the  model  in  which  many  small  particles  decelerate  and  distribute  themselves  throughout  a volume  of 
shocked  gas. 

A jet  mass  flux  of  25  to  50  percent  of  the  theoretical  maximum  seems  reasonable  in  view  of 
the  above  considerations.  This  would  lead  to  an  injection  of  15  to  30  grams  of  small  steel  particles 
into  the  driver  gas  for  each  driver. 

The  drivers  used  by  PI  in  Tests  1 to  4 do  not  exhibit  the  severe  jetting  characteristics  as 
do  the  drivers  used  in  Tests  5 to  8.  Results  of  Test  1 are  presented  in  Figure  22.  The  layer  of 
nitromethane  used  to  collapse  the  pressure  tube  was  increased  after  Test  4,  resulting  in  the  badly 
jetting  driver  operation  described  above.  It  is  not  certain  what  effect  the  limited  jetting 
observed  in  Tests  1 to  4 would  have  on  launcher  performance,  however,  with  the  severely  jetting 


Figure  21.  Schematic  diagram  of  jet  particle  trajectories  in  a jetting 
driver. 
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drivers  used  in  Tests  5 to  8,  launcher  performance  can  be  expected  to  be  adversely  affected.  The 
many  small  particles  may  cause  erosion  of  the  launcher  parts  downstream  of  the  injection  ports  and 
possibly  into  the  launch  tube.  If  driver  jet  particles  are  carried  into  the  launch  tube  they  can 
also  absorb  reservoir  energy  as  they  reaccelerate  with  the  flow. 

4.2  INJECTION  OF  SHOCKED  GAS  INTO  THE  COMPRESSOR  SECTION 

In  the  launcher  configuration  used  by  PI,  two  relatively  large  drivers  were  used  to  inject 
shocked  gas  into  the  compressor.  Because  the  combined  cross  sectional  area  of  the  drivers  is  some- 
what greater  than  that  of  the  compressor,  the  flow  of  injected  gases  can  choke.  The  resulting  re- 
duction in  mass  flowrate  can  be  a major  design  problem.  In  calculating  injection  gas  dynamics,  the 
most  accurate  t,,ation  of  state  for  hydrogen  should  be  used,  since  the  nature  of  the  injection  wave 
system  is  sensitive  to  flow  Mach  number,  isentropic  index,  etc.  In  addition,  a considerable  fraction 
of  the  driver  gases  may  be  trapped  between  the  driver  termination  flange  and  the  injection  ports  in 

the  PI  design,  even  if  there  is  no  reduction  in  injection  mass  flowrate. 

4.2.1  I ejection  Through  a Standing  Oblique  Shock 

In  the  geometries  tested  by  PI,  the  shocked  driver  gases  must  negotiate  a turn  during  injec- 
tion. In  Tests  1 through  7 the  turning  angle  was  45°  for  both  drivers.  In  Test  8,  only  one  driver 

was  used  and  the  turning  angle  was  reduced  to  24°.  For  Tests  1 through  7,  the  compressor  cross 
section  available  to  the  flow  was  less  than  the  sum  of  the  cross  sections  of  the  two  drivers. 
Specifically, 


A 

compressor  . q 
*both  drivers 

In  Test  8 the  cross  section  available  to  the  flow  increased  from  driver  to  compressor.  The  area 
ratio  was  1 .479. 

It  is  postulated  that  under  favorable  conditions  and  after  a start-up  process,  most  -of  the 
driver  gases  will  be  injected  through  a standing  oblique  shock  at  the  injector  port  as  shown  below: 
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Here 


e = shock  angle 


6 = flow  deflection  angle 


The  jump  relations  across  the  oblique  shock  are  solved  iteratively,  using  the  full  equation  of  state 
for  hydrogen  given  by  Equations  (A-j ) to  (A^)  and  noting  the  invariance  across  the  shock  of  the  tan- 
gential velocity  component  (Reference  5). 

The  flow  deflection  angle  is  computed  using 


6 = tan"' 


tan  Q (1  - P2/P3) 
Pg/p^  tan^O  + 1 


(10) 


In  Figure  23,  the  variation  of  deflection  angle  with  shock  angle  is  shown  for  the  range  of  shock  ve- 
locities 7.5  km/sec  < s < 9 km/sec.  The  maximum  possible  flow  deflection  angle  with  this  postulated 
wave  system  is  about  13°  with  a shock  velocity  of  7.8  km/sec,  characteristic  of  a properly  opera- 
ting driver.  If  covolume  effects  are  suppressed,  the  calculated  maximum  flow  deflection  angle  is 
21°,  which  illustrates  the  importance  of  accurately  representing  the  equation  of  state. 

By  the  above  analysis,  the  minimum  required  area  ratio  to  pass  the  flow  is  formed  from 


^3  _ fsin  (e  - (5) 

A„  [_  sin  6 

max 


For  a properly  operating  driver 


= 0.857 


min 


(11) 


It  must  be  conclude(f  that  a single  stationary  oblique  shock  cannot  turn  the  flow  in  the  PI 
launcher  configuration,  since  the  injection  angle  is  too  large  and  the  available  compressor  cross 
section  is  too  small.  Most  likely,  a normal  shock  must  propagate  back  up  the  driver  injection  ports 
to  reduce  the  flow  to  subsonic  conditions  for  injection. 


4.2.2  Subsonic  Injection  of  Driver  Gases 

A simple  flow  model  is  presented  here  to  provide  a rough  estimate  of  mass  flowrate  reduction 
as  the  driver  gas  is  injected  into  the  compressor.  After  a transient  startup  process,  a return 
shock  is  assumed  to  propagate  back  up  the  injection  port  into  the  explosive  driver.  The  flow  behind 
the  return  shock  is  subsonic  and  is  assumed  to  choke  as  it  flows  into  the  smaller  available  cross 
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Figure  23.  Influence  of  driver  shock  velocity  o,i  oblique  shock 
turning  angle  at  injection. 
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section  of  the  compressor.  The  strength  of  the  return  she:''  is  adjusted  so  that  the  flow  down  the 
compressor  is  at  sonic  or  choked  conditions  as  shown  below; 


This  simple  model  neglects  the  change  in  momentum  of  the  flow  as  it  turns  and  other  three- 
dimensional  effects,  such  as  a vena  contracta,  which  would  tend  to  further  restrict  mass  flowrates. 

Conditions  across  the  return  shock  are  calculated  iteratively  using  the  full  equation  of  state 
for  hydrogen.  The  subsonic  flow  behind  the  return  shock  (State  3a)  then  undergoes  an  area  reduction 
and  accelerates  to  sonic  conditions  (State  3b)  according  to  the  steady  energy  and  continuity  equa- 
tions and  the  isentropic  condition,  AS  = 0. 

For  example,  in  the  geometry  jf  PI  Tests  1 to  7,  the  conditions  behind  the  driver  shock  for 
ideal  driver  behavior  are: 


P2  = 0.974  kb 

U2  = 6-25  km/sec 
p2  = 0.0098  gm/cm’ 

Behind  a return  shock  propagating  back  up  the  driver  at  0.75  km/sec,  the  conditions  are: 

P3a  = 3.5  kb 
u^j  = 2.52  km/sec 
pjj  = 0.0210  gm/cm’ 

After  injection  into  the  compressor  section,  the  flow  is  sonic  with: 

P3b  = 2.1  kb 

U3b  =4.7  km/sec 
p3b  = 0.0153  gm/cm’ 
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This  combination  of  conditions  is  sufficient  to  provide  sonic  flow  in  the  cross  sectional  area  which 
is  available  in  the  compressor.  For  the  return  shock  speed  of  0.75  km/sec,  the  mass  flowrate  is  re- 
duced to  86  percent  of  its  unimpeded  value. 

The  compressor  cross  sectional  area  ratio  required  to  choke  the  flow  is  shown  in  Figure  24 
as  a function  of  return  shock  velocity.  The  reduction  in  mass  flowrate  as  a function  of  return 
shock  velocity  is  also  given.  From  this  figure,  it  can  be  shown  that  a vena  contracta  of  10  percent 
of  the  compressor  diameter  would  further  reduce  mass  flowrate  to  about  75  percent  of  its  unimpeded 
value.  For  precisely  timed  compressor  piston  arrivals  typical  of  the  PI  tests,  these  reduced  mass 
flowrates  could  result  in  significant  gas  losses  during  injection. 

4.2.3  Driver  Termination  Dynamics 

In  the  PI  launcher  configurations,  the  termination  of  the  explosive  drivers  is  considerably 
upstream  of  the  injection  ports.  Simple  calculations  indicate  that  the  rarefaction  generated  by 
the  termination  of  the  explosively  formed  piston  can  bring  to  rest  (and  therefore  trap)  a substan- 
tial fraction  of  driver  gases  within  the  drivers.  For  conditions  typical  of  PI  Test  7,  rough  esti- 
mates indicate  that  35  to  40  percent  of  the  driver  gas  can  be  lost  to  the  system  by  this  mechsnism 
(after  accounting  fur  the  return  shock  condition  considered  above). 

In  conclusion,  a substantial  fraction  of  driver  gases  was  not  injected  in  the  PI  tests  because 
of  explosive  driver  termination  dynamics  and  mass  flowrate  reduction  during  injection.  Thus,  the  ef- 
fective G/Mp  ratios  in  the  PI  tests  might  have  been  only  about  half  as  great  as  expected.  Since  PI 
launcher  performance  was  already  shown  to  be  sensitive  to  changes  in  the  G/Mp  ratio  (Section  3.4), 
it  is  not  surprising  that  muzzle  velocities  were  substantially  lower  than  anticipated. 

4.3  PROJECTILE  FRICTION 

Friction  between  the  projectile  and  launch  tube  can  substantially  reduce  projectile  accelera- 
tion, if  the  prelaunch  clearances  are  too  small.  A 5-kb  pressure  applied  to  the  base  of  an  aluminum 
projectile*  which  is  free  to  expand  in  the  radial  direction  dictates  a minimum  acceptable  prelaunch 
clearance  of  about  0.003  cm.  The  PI  clearances,  before  the  barrel  was  enlarged  due  to  erosion, 
ranged  from  approximately  0 to  0.005  cm.  Projectile  friction  clearly  existed  in  some  instances. 

For  example,  assuming  negligible  prelaunch  clearance,  the  radial  wall  stress  is  (Reference  4) 


1 - u. 


£ 

2 


(12) 


Modulus  of  elasticity  = 700  kb  (10’  psi),  Poisson's  ratio  : 0.33. 


where  = Poisson's  ratio 

p = projectile  base  pressure 

For  a projectile  of  diameter,  D,  and  length,  S,,  the  accelerating  force  is  p(Ti/4)/d^,  and  the 
retarding  force  is  f^[Mp/(l  - Mp)]  P/2,  where  f^  is  the  coefficient  of  sliding  friction.  The  ratio 
of  retarding  force  to  accelerating  force  is 


(13) 


This  ratio  is  independent  of  pressure  level.  For  the  PI  configuration,  the  projectile  H/d  is  2.4, 
tile  coefficient  of  sliding  friction  (aluminum  on  steel)  is  0.47,  and  Poisson's  ratio  is  0.33.  The 
ratio  of  retarding  force  to  accelerating  force  is  calculated  to  be  slightly  greater  than  1.  This 
is  clearly  not  possible,  but  this  simple  approach  does  indicate  that  projectile  friction  can  be  a 
significant  loss  mechanism  if  very  close  tolerances  are  used. 

4.4  BARREL  BOUNDARY  LAYER  EFFECTS 

The  shape  of  the  boundary  layer  in  a launcher  is  approximated  below. 


The  boundary  layer  has  a number  of  competing  effects  on  ballistic  performance  which  warrant 
the  following  general  discussion  before  proceeding  with  specific  evaluations. 

The  condition  at  the  wall  of  zero  velocity  causes  a retarding  effect  to  exist  through  the 
shear  action  within  the  boundary  layer.  The  retardation  approaches  zero  as  the  edge  of  the  boundary 
layer,  or  the  core  flow,  is  approached.  If  the  internal  surface  of  the  launcher  barrel  were  at  the 
same  temperature  as  the  core,  this  retardation  would  cause  a deficit  of  mass  flux  compared  to  the 
core  flow.  However,  local  continuity  considerations  would  dictate  that  the  core  be  accelerated  to 
compensate  for  the  deficit  of  mass  flux  within  the  boundary  layer.* 

The  internal  surface  temperature  of  the  launcher  barrel  is  usually  much  lower  than  the  core 
temperature.  Consequently,  the  higher  densities  near  the  wall  within  the  thermal  boundary  layer 
* 

It  is  more  precise  to  say  that  the  core  flow  would  behave  as  though  the  flow  were  inviscid,  but 
with  the  core  flow  area  varying  with  distance  down  the  barrel.  This  influences  global  continuity, 
which  may  or  may  not  cause  the  local  core  flow  to  be  accelerated. 


54 


tend  to  compensate  the  mass  flux  for  the  lower  velocities  within  the  hydrodynamic  boundary  layer. 
The  net  effect  is  that  the  core  flow  is  typically  influenced  very  little  by  the  boundary  layer. 
This  effect  depends  upon  the  magnitude  of  the  boundary  layer  displacement  thickness 


i 
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6*  = 


=i^i/ 


dy 


(14) 


where  p - density 
u = velocity 

y = direction  normal  to  surface 
& = boundary  layer  thickness 
i = core  property 

which  is  typically  near  zero  in  a launcher  barrel. 

If  the  core  flow  is  not  affected  by  the  boundary  layer  (i.e.,  6*  « 0),  then  the  ballistic 
performance  will  not  be  modified.  That  is,  in  spite  of  significant  heat  and  momentum  transfer  from 
the  gas  flow  behind  the  projectile,  the  influence  on  the  projectile  base  pressures  will  be  negligible 
as  long  as  the  boundary  layer  does  not  extend  to  the  center  of  the  barrel.  In  other  words,  boundary 
layer  momentum  and  heat  transfer  effects  could  be  ignored  in  terms  of  ballistic  performance,  if  the 
launcher  L/D  is  small  enough  to  obviate  boundary  layer  closure. 

Subsequent  to  boundary  layer  closure,  shear  is  transmitted  across  the  entire  cross  section 
(leading  to  reduced  base  pressures),  and  the  entire  cross  section  obtains  lower  temperatures  due  to 
heat  transfer  (also  leading  to  lower  base  pressures).  Consequently , the  discriminating  boundary 
layer  parameter  is  its  thiclkness,  6,  which  is  evaluated  below. 

4.4.1  Momentum  Losses  Due  to  Boundary  Layer  Closure 

The  nonsteady  boundary  layer  is  evaluated  herein  based  on  the  integral  analysis  of  Reference 

7.  The  details  of  the  evaluations  that  were  performed  for  this  program  are  presented  in  Reference 

8.  The  following  presentation  outlines  this  detailed  evaluation. 

The  nonsteady  integral  boundary  layer  momentum  equation  is  derived  for  quite  general  boundary 
conditions  in  Reference  7.  A restricted  form  of  this  equation  for  stationary  wall  conditions  and 
zero  transpiration  {or  ablation)  is; 
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i£  + 2 ^ ‘il  I F - -w 

at  G 3x  (G  at  G U.  dt  G ax'  > G' 


F = p^u.  e 


® = /o^  (1  - U/U.)  dy 

1 

H = 6*/9 


G = H - 6*/0 
P 

" fp^  ' p/Pi^  '^y 

= wall  shear  5 p^.u?  C^/2 
X = axial  coordinate 
t = time 

Reynolds  numbers  within  launcher  barrels  are  typically  high  enough  to  produce  turbulent 
boundary  layer  conditions  along  all  surfaces.  Linder  these  conditions,  but  for  steady  flow,  a 
convenient  relation  for  friction  factor  is  (Reference  9) 


V-  = 0.013  I 
2 P 


I = 

P P^  Vp2 


' symbol  denotes  properties  evaluated  at  reference  enthalpy  conditions,  h'  = 0.36  h^  + 0.19  h^  + 
0.45  h 

w 

h^  = recovery  enthalpy  = h.j  + R u^/2 

R = recovery  factor  (usually  slightly  less  than  1.0  in  gases) 
w denotes  gas  properties  evaluated  at  the  wall  temperature 
p = viscosity  (=  1.137  x 10'^  t‘"(°K)  - gm/(cm-sec) , u = 0.75) 


Various  boundary  layer  shape  parameters  have  been  evaluated  (Reference  10)  based  on  the  ap- 
proach outlined  in  Reference  9.  The  results  for  a turbulent  1/7  power  velocity  profile  are  pre- 
sented in  Figure  25.  It  is  noted  that  G is  approximately  equal  to  1.3,  independent  of  wall  temper- 
ature. Note  also  that  the  function,  1^,  increases  with  decreasing  wall  to  edge  temperature  ratio 
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(Ip  1.5  at  - 0.2).  Consequently,  the  shear  function  employed  in  the  remainder  of  this  evalu- 

ation is  based  on  Equation  (16)  and  these  latter  observations. 


! 0.013  I 

w _p 

_^i^i 

0.015  p. 

G ■ G 

piU.e^^ 

1 ^ ) 

1 h ) 

(17) 


Equation  (15)  can  be  solved  using  the  method  of  characteristics.  In  this  method,  a partial 
differential  equation  in  two  independent  variables  is  solved  by  simultaneous  solutions  of  two  total 
differential  equations  in  one  independent  variable  (one  equation  being  the  integral  momentum  equa- 
tion, the  other  being  the  equation  for  the  characteristic).  Choosing  time  as  the  independent  vari- 
able, after  inclusion  of  Equation  (17)  the  appropriate  equations  are: 


) H - G j ^ 

■(1  + Hj  H - G dX|  ■ 

3u. : 
1 ■ 

1 dt\ 

G G u^  dtj^ 

3x  ! 

0.015  p 


'A 


p . u . F 
1 ,c  1 ,c 
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(IB) 


dt  I 


c 


(19) 


where  c denotes  properties  along  a characteristic. 

Equations  (18)  and  (19)  can  be  solved  for  arbitrary  core  flow  condition  histories.  However, 
for  this  effort,  particular  core  behaviors  have  been  selected  for  certain  idealized  ballistic  cycles, 
such  that 


3^ 

3x 


c 


dt  ’ ^ 


C(t)  (to  be  defined) 


(20) 


The  solution  of  Equation  (IS)  can  be  obtained  by  formal  means  after  invoking: 

• Equation  (20) 

• Ihe  boundary  condition  of  zero  boundary  layer  thickness  at  the  origin  of  a characteris- 
tic (at  time  t.| ) 

• The  simplifications  that  H and  G are  constants. 


\ 5 '^1  '^2 
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(21) 
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Note  that  all  the  parameters  in  Equation  (21)  vary  with  time  along  characteristics.  The  character- 
istic is  developed  through  solution  of  Equation  (19).  The  properties  and  the  solution  of  Equation 
(19)  depend  upon  the  core  flow  descriptions  which  are  evaluated  in  the  following  paragraphs. 


I 
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Two  limiting  ballistic  cycles  are  considered; 
t Constant  base  pressure  launcher  (CBPL) 

• Simple  wave  launcher  (SWL) 

It  is  convenient  to  work  in  nondimensional  coordinates  to  describe  ballistic  behavior 


X* 


2^0 


t*  = t 

2*. 


II*  = 


_L)_ 


2a 


7 ; isentropic  exponent 
V ; initial  acceleration  of  the  projectile 

a = initial  speed  of  sound  in  the  chamber  i 

0 

Equation  (21)  is  transformed  with  these  variables  to:  j 


where  6*  ^ ^ 9 

2^0 

0 denotes  conditions  at  the  initiation  of  projectile  motion 

z denotes  the  origin  of  a characteristic,  either  at  the  chambrage  (z  = c),  or  along  the 
barrel  where  the  characteristic  is  originated  at  the  time  of  projectile  passage  (z  = p) 
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! 
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In  the  constant  base  pressure  launcher,  acceleration  is  always  equal  to  the  initial  accelera- 
tion, and  so  the  projectile  motion  is  easily  evaluated.  Continuity,  inomentum,  and  energy  considera- 
tions dictate  that  gas  velocity  be  independent  of  distance  along  the  barrel,  but  vary  with  time  in 
the  same  way  as  the  projectile  velocity  does.  (This  solution  is  presented  by  Seigel  in  Reference  2.) 


uT  = t* 


(23) 


This  cycle  is  particularly  simple  in  that  characteirstics  vary  in  the  same  way  throughout  the  barrel, 
and  the  function,  C.  is  invariant  (C  - 1.0). 

A simple  wave  launcher  with  the  chamber  and  barrel  cross  sections  equal  is  analyzed.  The  lo- 
cation of  the  chamber  head  end  is  chosen  to  be  far  enough  away  from  the  chambrage  plane  so  that  the 
reflection  of  the  lead  expansion  wave  (which  is  generated  at  the  beginning  of  projectile  motion) 
does  not  cross  the  chambrage  plane  prior  to  completion  of  the  ballistic  event. 

This  cycle  is  also  considered  by  Seigel  (Reference  2),  but  without  obtaining  an  explicit 
evaluation  of  the  distributions  of  gas  properties  along  the  barrel.  Reference  8 has  shown  these 
properties  to  vary  isentropically,  with  the  following  velocity  distribution. 


X*  + t* 
1 + t* 


(24) 


such  that 


_2 

CgwL  = n + (25) 

The  equations  of  projectile  motion  and  all  equations  which  are  necessary  to  the  evaluation  of 
Equation  (22)  are  presented  in  Table  5.  Both  cycles  yield  pressures  which  decrease  with  distance 
down  the  barrel.  However,  pressure  decreases  with  time  at  the  chambrage  plane  of  the  SWL,  whereas 
it  must  increase  with  time  to  maintain  base  pressure  in  the  CBPL.  Gas  temperatu'"es  change  in  the 
same  directions  as  the  pressures,  requiring  isentropic  compression  of  the  chamber  gases  with  time 
for  the  CBPL. 

In  terms  of  two-stage  launchers,  the  CBPL  can  be  thought  of  as  representing  a specific  his- 
tory of  compressor  piston  movement  toward  the  chambrage  plane  during  the  ballistic  event.  The  SWL, 
however,  can  be  thought  of  as  representing  a specific  piston  movement  away  from  the  chambrage  plane. 
Presumably,  results  evaluated  from  these  idealized  cycles  will  reveal  representative  limits  of 
boundary  layer  behaviors  for  two-stage  launchers. 
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TABLE  5.  TIME  (t*)  VARIATIONS  OF  PROJECTILE  MOTION  AND  INVISCID  CHARACTERISTIC 
PROPERTIES 


Parameters 


(1  + t*) 


1 0 - -T.z  ^ 


(1),  (2),  (3):  See  notes  on  following  page 
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NOTES  TO  TABLE  5 


(1)  t*  > 


p denotes  projectile 

c denotes  chambrage  plane 

z denotes  that  either  p or  c applies 

(2)  K*  = X*  (t*  ),  projectile  originated  characteristic 

1 1 p p I » P 

(3)  X*  = 0,  chambrage  plane  originated  characteristic 

I } c 
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Projectile  motions  are  presented  for  both  cycles  in  Figure  26.  These  are  compared  with 


r 


muzzle  parameters  which  have  been  evaluated  in  Section  3.*  For  equivalent  muzzle  velocities  and 
initial  accelerations,  the  SWL  travel  time  and  barrel  length  are  considerably  longer  than  for  the 
CBPL.  Howeve  for  the  same  muzzle  velocities  and  barrel  lengths,  the  travel  time  is  shorter  for 
the  SWL,  requiring  higher  initial  accelerations.  These  comparisons  are  made  for  a muzzle  velocity 
of  u*  = 1.0  in  the  following  example: 


COMPARISON  OF  CBPL  AND  SWL  BALLISTICS,  u*  = 1.0 

m 


t* 

ni 

L*(=  X*) 

/^WL  \ 

\^CBPL/u,u* 
' m 

\4bPL/v,u* 

/ ^_S_WL  \ 
Vcbpl/l,u* 

/ ^SWL  \ 
V-'CBPL/L,u* 

SWL 

1.99 

1.0 

1 .99 

2.44 

0.82 

2.44 

CBPL 

1.0 

0.41 

The  simplified  boundary  condition  of  zero  boundary  layer  thickness  at  the  chambrage  plane 
for  all  times  is  not  strictly  correct.  It  would  be  a reasonable  approximation  for  a design  with  a 
sharp  transition  between  the  chamber  and  the  barrel,  but  in  terms  of  heat  transfer  this  would  be 
poor  design  practice.  In  steady  flow,  boundary  layers  are  thinned  under  accelerating  inviscid  flow 
conditions.  Since  flow  would  be  accelerated  for  practical  transition  section  designs,  it  is  proba- 
ble that  at  least  the  boundary  layers  are  thinnest  at  the  chambrage  plane.  The  assumption  of  zero 
thickness  here  has  a negligible  effect  on  boundary  layer  closure  evaluations,  and  a conservative 
effect  in  terms  of  heat  transfer  evaluations. 

When  H is  set  equal  to  0.2  G,  B.|  becomes  -1.  In  this  case,  the  simplicity  of  the  CBPL  al- 
lows a formal  integration  to  be  carried  out  in  Equation  (22).  (The  integration  must  be  performed 
by  numerical  means  for  the  SWL  due  to  the  complexity  of  the  equation  for  the  characteristics.)  The 
final  relation  for  momentum  thicknesses  for  the  CBPL  is  (H  = 0.2  G). 


★ 

The  SWL  and  CBPL  results  generally  span  the  earlier  results.  The  results  of  Section  3 cluster 
near  the  SWL  resQlts,  implying  the  possibility  of  higher  ideal  muzzle  velocities  than  found  in 
Section  3. 
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Fiqure  26.  Generalized  ballistic  performance. 
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where 


A = 


R - 2.(A-_JJ 
^3  G 


(y  + 1)^ 


For  both  cycles,  maximum  momentum  thickness  occurs  very  near  the  characteristic  originated  at 
zero  time  which  is  common  to  both  the  chambrage  plane  and  projectile.  For  this  characteristic, 
Equation  (26)  becomes: 
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(1  + Bjt*  ) - 1 


(1  + B3t*  ) 


1 

Y-1 


(27) 


For  a very  high  speed  of  sound  in  the  chamber,  t*  is  always  small,  leading  to  the  following  approxi- 
mation of  Equation  (27): 


V V 


0*  = 0.0315  I 

U(Y^l)Xy 


0 \ 


(28) 
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Example  solutions  of  Equation  (22)  along  the  common  characteristics  have  been  performed  using  both 
ideal  cycles  and  a muzzle  velocity  of  6 km/sec.  The  barrel  length  is  constrained  for  these  calcula- 
tions to  300  diameters,  assuming  a 2.54-cm  diameter  barrel.  For  the  sake  of  consistency  with  the 
derivation  of  Equations  (26)  and  (27),  a wall  temperature  ratio  of  0.55  has  been  selected  which 
yields  H = 0.2  G (Figure  25),  and  a value  of  6/6  = 8.6.  For  these  assumptions,  closure  occurs  when 
9 is  equal  to  0.148  centimeters. 

The  results  are  presented  in  Figure  27.  The  boundary  layer  is  thicker  for  a given  projectile 
position  for  the  CBPL  than  it  is  for  the  SWL.  This  follows  from  the  constraint  of  the  same  muzzle 
velocity  and  barrel  length  for  both  cycles.  The  SWL  starts  with  a higher  chamber  pressure  and  ob- 
tains higher  projectile  velocities  at  all  positions  along  the  barrel  except  at  the  muzzle.  Gas  tem- 
peratures are  lower  in  the  SWL,  so  viscosities  are  lower.  For  the  early  part  of  the  SWL  nrotion,  the 
gas  densities  are  higher  than  in  the  CBPL.  The  values  of  these  parameters  are  such  that  core  Rey- 
nolds numbers  at  the  point  of  peak  boundary  layer  thickness  are  always  higher  for  the  SWL  than  the 
CBPL  (i.e.,  under  the  constraints  of  same  muzzle  velocity  and  barrel  length).  Consequently,  the  SWL 
boundary  layers  are  thinner. 
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Parametric  values  of  H have  been  employed  in  the  calculations  for  the  SWL.  Although  peak 
momentum  thicknesses  decrease  with  tl.  6/e  increases  with  H.  Consequently,  the  location  of  the  pro- 
jectile when  closure  is  obtained  is  almost  independei.t  of  H (projectile  location  varies  less  than  15 
centimeters  in  Figure  27  for  the  range  of  H considered). 

Evaluations  of  the  exact  solution  of  Equation  (22),  (i.e..  Equation  (26)  with  H = 0.2  G)  and 
the  approximate  solution  (Equation  '28)),  for  the  CBPL  are  presented  in  Figure  27.  These  results 
compare  within  about  3 percent  near  the  point  of  closure,  so  the  approximate  results  are  sufficiently 
accurate  for  the  following  purposes. 

It  is  of  interest  to  determine  the  effects  of  muzzle  velocity,  chamber  temperature,  and 
launcher  barrel  L/D  or  the  location  of  the  projectile  at  the  time  of  closure.  This  can  be  approxi- 
mated from  Equation  (28)  through  retting  = 2(5/11)  ’ 


28.7 

^ 

V .a'  m 

\ 0 P > 


where  m^  = projectile  mass  per  unit  area. 

Using  fonner  dependencies  of  viscosity  and  sound  speed  on  temperature,  and  a perfect- ideal  gas, 


(»c. 


osure  “ VD 


Equation  (29)  has  been  evaluated  for  the  projectile  areal  density  that  was  used  in  the  PI 
launcher  (17  gm/cm^),  and  a 6/0  of  7.1.  The  results  are  presented  in  Figure  23. 

It  is  significant  to  note  that  if  insufficient  mass  is  injected  into  the  chamber,  the  ideal 
muzzle  velocity  will  be  low.  Figure  20  shows  that  when  this  happens,  boundary  layer  closure  occurs 
at  smaller  values  of  projectile  travel.  This  in  turn  makes  the  attainment  of  the  ideal  muzzle  ve- 
locity difficult,  since  performance  falls  off  after  boundary  layer  closure.  Thus,  it  is  not  only 
important  to  realize  full  chamber  conditions  in  terms  of  obtaining  maximum  ideal  velocity,  but  it 
is  also  important  in  order  to  minimize  deoartures  from  ideal  behavior  due  to  boundary  layer  effects. 

The  results  for  the  SWL  in  Figure  28  show  a greater  sensitivity  of  closure  to  ideal  muzzle 
velocity  than  the  CBPL  does.  It  is  probable  that  the  SWL  results  approach  the  CBPL  results  asymptoti- 
cally with  decreasing  muzzle  velocity,  based  on  the  observed  approach  of  the  SWL  ballistic  performance 
results  to  the  CBPL  results  at  smaller  u*  (Figure  26). 
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Although  the  SWL  allows  more  projectile  motion  before  closure,  it  is  important  to  remember 
that  this  is  being  obtained  here  by  virtue  of  higher  initial  accelerations.  If  both  guns  are  con- 
strained to  the  same  initial  acceleration,  the  higher  L/D  required  for  the  SWL  offsets  this  apparent 
advantage  (applying  the  CBPL  rules  of  Equation  (28)  to  the  SWL,  doubling  the  gun  L/D  will  probably 
reduce  the  projectile  travel  to  closure  for  the  SWL  by  about  26  percent).  The  increased  L/D  de- 
creases (x  /D)  , , and  increases  the  distance  that  must  be  travelled  after  closure  occurs. 

' p closure 

The  results  of  these  evaluations  can  now  be  summarized: 

• Boundary  layers  do  not  significantly  influence  ballistic  performance  until  boundary 
layer  closure  occurs 

• Launcher  L/D  should  be  small  to  minimize  boundary  layer  effects  on  ballistics 

t Allowable  L/D's  increase  with  increasing  ideal  muzzle  velocity  (i.e.,  increasing  accelera- 
tions) 

• Phenomena  which  degrade  the  chamber  conditions  in  such  a way  as  to  reduce  the  ideal 
muzzle  velocity  also  reduce  the  opportunity  to  achieve  even  the  reduced  ideal  muzzle 
velocity  because  of  adverse  boundary  layer  effects 

• Quantitative  and  accurate  evaluations  of  the  effects  of  boundary  layer  closure  on  per- 
formance require  that  momentum  and  heat  transfer  be  coupled  to  the  computation  of  bulk 
behavior 

• The  barrel  L/D  employed  in  the  PI  launcher  tests  was  much  greater  than  is  useful.  How- 
ever, the  excess  barrel  length  probably  neither  enhanced  nor  reduced  the  muzzle  velocity 

Estimates  of  the  effects  of  heat  transfer  on  ballistic  performance  are  presented  in  the  next 
subsection. 

4.4.2  Energy  Losses  Due  to  Heat  Transfer 

It  is  desirable  to  estimate  the  upper  limits  of  ballistic  losses  due  to  heat  transfer.  In 
addition,  heat  transfer  must  be  evaluated  to  estimate  the  thermal  response  of  the  launcher  barrel 
internal  surface  material.  Consequently,  the  analysis  in  Section  4.4.1  is  extended  here  and  will 
consider  the  spatial  and  temporal  variations  of  heat  transfer,  as  well  as  the  effects  of  heat  trans- 
fer on  ballistic  performance. 

It  is  common,  and  appropriate  in  approximate  evaluations  of  this  sort,  to  invoke  Reynolds' 
analogy  to  evaluate  heat  transfer,  that  is; 


69 


Operating  on  Equation  (31)  using  the  transformed  variables,  a heat  transfer  coefficient  can 
be  evaluated  as  follows: 


Example  axial  distributions  of  heat  transfer  coefficient  obtained  through  evaluations  of 
Equation  (32)  for  both  the  CBPL  and  SWL  cycles  are  presented  in  Figure  29(a).  The  times  chosen  for 
presentation  represent  times  when  the  projectile  would  reach  the  muzzle  for  the  conditions  noted  on 
the  figure.  The  infinite  heat  transfer  coefficients  at  both  the  base  of  the  projectile  and  at  the 
chambrage  plane  are  unrealistic.  The  finite  thickness  of  the  boundary  layer  at  the  chambrage  plane 
was  discussed  in  Section  4.4.1.  Corrections  to  this  result  can  by  obtained  by  theoretical  means 
quite  easily  in  principle,  but  the  effort  is  beyond  the  scope  of  the  present  program.  However,  except 
at  the  projectile,  the  results  in  Figure  29(a)  at  values  of  x/L  greater  than  about  0.05  are  probably 
quite  accurate.  As  will  be  seen  later,  the  condition  at  the  base  of  the  projectile  is  unimportant, 
because  the  attenuation  of  the  heat  flux  there  is  very  rapid. 

Minimum  heat  transfer  coefficients  exist  near  the  points  of  maximum  boundary  layer  thickness 
(x  ' 580  cm).  Transfer  coefficients  for  the  CBPL  are  on  the  order  of  5 to  30  percent  higher  than 
for  the  SWL.  Note  that  these  results  are  somewhat  in  error  (probably  high)  in  the  region  where  the 
boundary  layer  is  fully  developed  (between  about  150  and  750  cm  for  both  cycles).  Based  on  the 
higher  heat  transfer  coefficients  and  recovery  enthalpies*  (Figure  29(a)  and  29(b))  at  the  time  of 
projectile  exit,  the  CBPL  heat  transfer  rate  is  considerably  higher  than  that  for  the  SWL. 

The  critical  area  for  heat  transfer  is  near  the  chambrage  plane,  where  transfer  coefficients 
are  high  and  the  surface  is  exposed  to  hot  gases  for  the  longest  times.  Transfer  coefiicient  his- 
tories at  two  locations  near  the  chambrage  plane  for  both  cycles  are  presented  in  Figure  30(a).  The 
rapid  attenuation  of  transfer  coefficient  after  passage  of  the  projectile  and  the  attenuation  of 
heat  transfer  with  dist  ce  (Figure  29(a))  are  apparent  in  this  figure. 


Recovery  enthalpies  have  been  computed  employing  a recovery  factor  of  unity. 
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It  is  somewhat  surprising  that  the  average  transfer  coefficients  between  the  time  of  pro- 
jectile passage  and  the  tine  the  projectile  leaves  the  barrel  are  about  the  same  for  both  idealized 
cycles.  The  lengths  of  the  heating  pulses  are  also  approximately  the  same  between  cycles.  However, 
considering  the  time  histories  of  local  recovery  enthalpy  shown  in  Figure  30{b)  (the  differences  in 
recovery  enthalpy  between  the  CBPL  and  SWL  increase  with  muzzle  velocity),  important  differences  in 
heat  transfer  exist  late  in  the  propulsion  cycle. 

The  time  variations  of  heat  transfer  rates  {based  on  a constant  barrel  temperature  of  1,000°K*) 
over  the  internal  surfaces  (i.e.,  integral  of  heat  fluxes  over  the  exposed  surface)  are  shown  in 
Figure  31.  The  principal  cause  for  the  rate  increase  with  time  is  the  increase  in  exposed  surface 
area  with  time  rather  than  heat  flux. 

Note  that  these  predicted  heat  transfer  rates  are  conservatively  high  after  boundary  layer 
closure.  This  is  primarily  because  heat  loss  from  the  gas  causes  a reduction  of  heat  transfer  driv- 
ing potential  after  closure.  (The  curve  labeled  "Reference"  in  Figure  31  addresses  this  considera- 
tion, and  will  be  discussed  later.) 

Total  heat  loads  (time  integrals  of  rates  in  Figure  31)  for  representative  conditions  are 
presented  in  Figure  32.  These  loads  are  on  the  order  of  the  projectile  kinetic  energies  at  the 
muzzle  (kinetic  energy  = 3.7  x 10^  calories  for  Mp  = 86  gm,  u = 6 km/sec).  Consequently,  it  is 
clear  that  they  need  to  be  considered  in  terms  of  ballistic  performance. 

The  heat  load  in  Figure  32  for  the  CBPL  is  80  percent  larger  than  that  of  the  SWL  at  6 km/sec. 
The  heat  loads  for  the  two  cycles  approach  each  other  at  lower  velocities,  which  is  to  be  expected 
(based  on  equivalence  of  the  cycles  at  low  u*). 

Two  approaches  have  been  employed  to  estimate  the  effects  of  heat  transfer.  In  the  first 
approach,  total  heat  remov»l  ideally  takes  place  in  the  chamber  prior  to  projectile  motion.  Based 
on  early  results,  heat  losses  for  this  analysis  have  been  approximated  by  the  following  relationship; 

Q(cal)  = 3.02  X 10“  (u.^  (33) 

Results  of  evaluating  Equation  (33)  are  slightly  higher  than  those  for  the  SWL,  as  shown  in  Figure 
32. 

In  the  calculation,  the  reservoir  state  after  the  heat  loss  is 


^R,Q  ^R,Q=o  ■ ^ 

Obviously,  heat  losses  will  be  lower  for  higher  barrel  temperatures. 


(34) 
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and  the  modified  reservoir  parameters  are  computed  as  in  Section  3 for  dissociated  hydrogen,  includ- 
ing covolume  effects.  Subsequently,  to  obtain  the  estimated  performance  with  heat  transfer,  ideal 
performances  have  been  calculated  using  the  Seigel  curves.  The  results  are  presented  in  Table  6 
where  it  is  seen  that  the  maximum  effect  is  a reduction  of  less  than  9 percent  in  muzzle  velocity. 

In  an  attempt  to  distribute  the  heat  load  during  the  ballistic  cycle  in  a more  representative 
manner,  another  idealized  cycle  was  calculated  with  heat  transfer  directly  coupled  into  the  expansion 
process.  In  this  technique,  nonsteady  flow  effects  are  ignored  in  favor  of  assuming  uniform  spatial 
distribution  of  density,  pressure,  and  temperature.  In  addition,  convection  from  the  gas  is  evaluated 
by  applying  a fully  developed  turbulent  pipe  flow  heat  transfer  relationship.  The  analytical  tech- 
nique is  developed  in  Reference  11  for  application  to  a low  speed  pneumatic  propulsion  concept. 

Without  heat  transfer,  this  analysis  yields  higher  muzzle  velocities  than  that  obtained  using  the 
Seiqel  curves. 

The  condition  considered  for  application  of  this  technique  is  Case  No.  3 of  Table  6 (this  is 
also  Case  Mo.  7 of  Tables  3 and  4).  In  the  calculation,  system  volumes  and  masses  are  the  same  as 
in  Table  6.  However,  the  method  of  Reference  11  considers  chamber  and  barrel  cross  sections  to  be  the 
same  (this  requires  a very  long  chamber  for  the  present  calculation),  and  evaluates  heat  transfer  in 
the  chamber  as  well  as  in  the  barrel.  Because  of  the  low  velocities  in  the  actual  compressor  section, 
chamber  heat  losses  are  exnected  to  be  relatively  small  compared  to  losses  in  the  barrel.  This  consid- 
eration (in  addition  to  the  expected  minimum  ballistic  effects  of  heat  transfer  before  closure)  sug- 
gests that  the  estimated  heat  loss  effects  on  ballistics  will  be  conservatively  high. 

The  evaluation  is  based  on  a perfect-ideal  gas  which  has  been  represented  by  choosing  a suit- 
able fictitious  molecular  weight.  The  choice  was  based  on  matching  the  isentropic  exponent,  chamber 
density,  chamber  volume,  pressure,  and  temperature  of  Table  6. 

The  time  history  of  heat  transfer  rates  from  this  evaluation  is  presented  as  the  curve  labeled 
"Reference"  in  Figure  31.  The  rates  are  very  much  higher  than  the  former  data  in  Figure  31  during 
the  early  portion  of  the  cycle  due  to  the  computed  convection  within  the  chamber;  later,  the  rates 
decrease  due  to  the  loss  of  energy  from  the  gas.  The  history  of  the  heat  losses  for  the  CBPL  and  SWL 
are  considerably  more  accurate,  at  least  prior  to  the  times  of  boundary  layer  closure.  The  total 
heat  load  is  shown  in  Figure  32  for  the  computed  ideal  velocity  of  5.5  km/sec  (computed  from  the 
method  of  Reference  11).  This  heat  load  is  approximately  30  percent  higher  than  for  the  CBPL  at  the 
same  muzzle  velocity. 

The  result  in  Table  6 shows  an  ideal  muzzle  velocity  approximately  13  percent  lower  than  that 
computed  by  the  method,  shown  in  Figure  33.  The  computed  velocity  degradation  due  to  heat  transfer 
is  approximately  16  percent.  Also,  the  computed  pressure  and  temperature  at  the  time  of  projectile 
exit  are  significantly  influenced  by  heat  transfer,  as  shown  in  Table  7. 


TABLE  7.  ESTIMATED  EFFECTS  OF  HEAT  TRANSFER  ON  FINAL  GAS  PRESSURE 
AND  TEMPERATURE  (DATA  OF  FIGURE  33) 


If  the  calculation  is  extended  to  greater  barrel  lengths,  projectile  velocity  will  still  in- 
crease. However,  the  rate  of  increase  is  quite  diminished  at  an  L/D  of  300. 

Based  on  all  these  estimates,  it  is  concluded  that  the  heat  losses  in  the  PI  launcher  are  not 
a major  cause  for  performance  shortfall,  if  performance  is  judged  at  projectile  travels  of  less  than 
about  300  diameters  (a  shortfall  of  approximately  10  percent  is  to  be  expected  here).  If  the 
design  requires  additional  L/D  to  get  performance  (based  on  ideal  calculations),  the  expected  re- 
sults will  not  be  achieved.  In  the  PI  tests,  maximum  projectile  velocities  were  achieved  at  less 
than  430*  diameters  for  Test  0,  which  is  not  surprising.  Muzzle  velocity  in  Test  7 was  achieved 
in  about  200  diameters. 
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SECTION  5 


MATERIAL  RESPONSE  TO  HEAT  TRANSFER 

Heat  transfer  rates  to  the  barrel  are  extremely  high  Tor  the  muzzle  velocities,  projectile 
areal  densities,  and  the  gas  enthalpy  levels  of  interest  here.  Examples  of  spatial  and  temporal 
distributions  of  heat  transfer,  as  well  as  integrated  heating  loads  have  been  presented  in  Section  4.4. 
In  Section  5,  various  candidate  barrel  (or  barrel  insert)  material  responses  to  the  heating  encoun- 
tered during  specific  ballistic  cycles  are  discussed. 

The  evaluations  of  material  response  were  performed  using  the  one  dimensional  “Aerotherm 
Charring  Material  Thermal  Response  and  Ablation  Program"  (CMA).  This  code  is  described  in  Ref- 
erence 12.  Input  to  the  CMA  code  include  pressure  and  recovery  enthalpy  (approximated  by  the  total 
enthalpy)  histories  from  the  ballistic  cycle,  and  the  heat  transfer  coefficient  distribution  from 
the  boundary  layer  calculations.  This  code  uses  real  gas  thermochemistry  (in  this  case,  for  hydro- 
gen) generated  with  the  "Aerotherm  Chemical  Equilibrium"  (ACE,  Reference  13)  code,  to  determine  the 
wall  enthalpy,  h^,  as  a function  of  wall  temperature  and  pressure.  Thus,  the  heat  flux  which  is  a 
boundary  condition  on  the  conduction  problem  depends  upon  its  solution  through  the  wall  temperature. 

The  transient  internal  conduction  computation  uses  a numerical  finite  difference  procedure. 

The  results  of  Subsection  4.4.2  suggest  that  heat  losses  to  the  barrel  do  not  have  a large 
impact  on  ballistic  performance.  Accordingly,  the  heat  transfer  problem  in  the  launcher  is  not  a 
problem  of  containing  the  heat,  but  rather  one  of  handling  the  heat  transfer  with  minimum  damage  to 
launcher  components . 

The  heat  transfer  can  be  handled  in  at  least  three  ways  by: 

• Selecting  a barrel  or  barrel  liner  material  with  properties  such  that  it  does  not  reach 
failure  temperature  before  completion  of  the  ballistic  event.  This  approach  makes  use 
of  the  heat  sink  capabilities  of  the  material 

• Actively  cooing  the  material  below  failure  temperature  (e.g.,  transpiration  cooling) 

• Allowing  the  material  to  fail  (or  ablate),  and  replacing  it  as  required 
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The  best  approach  is  the  application  of  heat  sink  materials,  since  they  are  passive  and 
(ostensibly)  will  not  need  to  be  replaced.  Active  cooling  may  be  costly  and  complex,  and  if  used, 
should  be  designed  to  cool  to  just  below  the  material  fail  temperature  to  minimize  system  convective 
heat  losses.  Ablative  liners  would  recede  intentionally  in  a controlled  manner  and  perhaps  would 
survive  several  launch  cycles.  This  effort  has  concentrated  on  evaluating  the  feasibility  of  heat 
sink  materials;  the  only  failure  mechanism  considered  was  melting. 


The  relative  perfomiance  potential  for  candidate  materials  can  be  predicted  using  the  results 
of  transient  conduction  in  a constant  property  semi-infinite  slab  subjected  to  constant  heat  flux 
(see  for  example.  Reference  14).  For  this  case  the  time  to  melt  is 


t 


m 


(35) 


where  AT  - 


k = thermal  conductivity 
c = specific  heat 

Among  structural  metals,  tungsten  has  the  longest  time  to  melt  for  a given  heat  flux,  be- 
cause of  its  extremely  high  melting  point.  On  the  other  hand,  steel  is  inexpensive  and  strong,  and 
copper  has  a high  thermal  conductivity.  Relative  to  tungsten,  these  materials  have  melting  times 
(at  constant  heat  flux)  of 

t = -L  t 

m,  steel  15  m,w 

and  t = T t 

m,  copper  4 m,w 

In  the  present  case,  an  additional  benefit  can  be  expected  from  high  melt  temperatures  through  the 
reduction  of  the  heat  flux  potential,  h^  - h^.  The  results  presented  below  demonstrate  that  this 
reduction  is  an  important  consideration. 

It  has  been  found  that  even  the  least  severe  heating  environment  can  be  expected  to  melt 
ordinary  heat  sink  materials.  Material  comparisons  were  made  from  thermal  response  predictions 
computed  for  the  3 km/sec  CBPL  heat  transfer  coefficient  and  recovery  enthalpy  histories  shown  in 
Figure  30,  together  with  the  appropriate  pressure  histories.  These  results  are  shown  in  Figure  34. 
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Figure  34.  Thermal  response  of  three  barrel  surface  materials  to  convection 
near  the  chambrage  plane.  CBPL,  conditions  of  Figure  30  (i.e., 

3 km/sec,  p = 1000  atm). 
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Relative  thermal  performances  (Figure  34)  are  approximately  as  expected.  These  predictions  are  sup- 
ported by  the  fact  that  melting  of  steel  occured  in  the  PI  launcher  at  the  indicated  muzzle  velocity.* 
Steel  and  copper  were  judged  to  be  inadequate  heat  sink  materials  for  the  desired  application 
(6  km/sec),  and  predictions  at  high  muzzle  velocities  were  restricted  to  tungsten. 

Tungsten  thermal  response  predictions  were  performed  for  a 6 km/sec  muzzle  velocity  at  two 
locations  in  a 2.54  cm  diameter  launcher  operating  either  as  a CBPL  or  as  a SWL.  The  heat  transfer 
coefficient  histories  are  similar  in  shape  to  those  shown  in  Figure  30,  but  have  a magnitude  5 or 
6 times  as  large.  The  thermal  response  predictions,  which  are  displayed  in  Figure  35,  indicate  that 
tungsten  will  melt  in  a CBPL  at  these  conditions,  but  will  not  melt  in  a SWL.  The  reason  for  the 
vast  difference  between  the  predictions  for  the  two  cycles  is  apparent  from  the  total  enthalpy  dis- 
tributions shown  in  Figure  36. 

The  total  enthalpy  is  defined  (for  a perfect  gas)  as 
h = Cp  T + U V2 

Cp  = specific  heat  at  constant  pressure 

At  the  time  the  projectile  leaves  the  muzzle,  the  velocity  at  the  chambrage  plane  is  6 km/sec  in  the 
CBPL  and  only  3.5  km/sec  in  the  SWL.  Also,  the  local  recovery  temperature  increases  in  a CBPL  and  de- 
creases in  a SWL.  The  net  result  is  that,  at  muzzle  time,  the  ratio  of  total  enthalpies  for  the  two 
cycles  is  about  1.8.  Since  the  enthalpy  of  hydrogen  at  the  melting  temperature  of  tungsten  is  almost 
as  large  as  the  SWL  total  enthalpy,  the  differences  in  total  enthalpy  between  the  two  cycles  has  a 
drastic  effect  on  heat  transfer.  This  effect  would  be  decreased  for  higher  initial  chamber  tempera- 
ture and  magnified  for  lower  temperatures. 

Although  the  SWL  results  appear  encouraging  from  a material  protection  point  of  view,  it  must 
be  kept  in  mind  that  base  pressures  are  prohibitively  high  at  early  times  in  this  mode  of  operation. 
Figure  37  describes  the  pressure  histories  at  the  chambrage  plane  (pp)  and  on  the  projectile  base 
(Pp)  required  to  achieve  the  projectile  muzzle  velocity  of  6 km/sec  for  this  example.  The  CBPL  pro- 
vides the  best  ballistic  performance  for  a given  base  pressure,  but  has  prohibitively  high  heat 
transfer  rates.  Thus,  it  is  apparent  that  alterations  to  a given  cycle  which  tend  to  solve  one  of 


* 


No  calculations  were  performed  for  this  case  for  a SWL.  Since  SWL  heat  transfer  rates  are  lower 
near  the  end  of  the  cycle,  it  is  probable  that  both  copper  and  steel  would  survive  for  this  ex- 
ample at  X = 30  inches  (76.2  cm),  but  neither  would  survive  at  x = 3 inches  (7.62  cm). 
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Figure  36.  Stagnation  enthalpy  history  input  to  calculation  of  Figure  35. 
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Figure  37.  Pressure  histories  associated  with  an  86  gram  projectile  muzzle  velocity 
of  6 km/sec,  762  cm  barrel  length. 
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these  problems  will  only  tend  to  aggravate  the  other.  However,  a compromise  may  be  feasible,  de- 
pending on  the  degree  of  control  of  compressor  piston  motion.  The  desired  motion  will  achieve  near 
constant  base  pressure  until  midway  in  cycle,  and  reduced  base  pressure  thereafter  in  the  manner 
of  a SWL  to  limit  recovery  enthalpies  and  heat  transfer  coefficients. 

It  may  be  that  a compromise  between  high  muzzle  velocity  at  limited  base  pressure  and  limited 
heat  transfer  cannot  be  found  which  will  allow  heat  sink  absorption  of  the  energy  without  melting. 
Thermal  protection  using  a subliming  ablator  (for  example,  carbon)  may  be  feasible  and  should  be 
carefully  considered  in  any  future  design  of  high  performance  light  gas  launchers.  The  necessary 
evaluation  is  beyond  the  scope  of  the  present  work. 
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SECTION  6 


RECOMMENDATIONS 

Aerotherm  recommends  that  the  hybrid  launcher  be  redesigned  to  achieve  the  desired  ballistic 
performance.  The  redesign  must  address  the  following: 

• Hydrogen  is  the  best  choice  for  the  working  medium.  Within  the  operating  range  of  the 
hybrid  launcher  (1  to  10  kilobars  and  1,000  to  6,000  degrees,  Kelvin),  it  is  imperative  to 
account  for  molecular  covolume,  vibrational  excitation,  and  dissociation. 

• Explosive  drivers  should  be  designed  to  operate  without  jetting  or  leaking 

t The  configuration  must  be  such  that  complete  driver  gas  injection  is  realized.  It  may 

be  necessary  to  use  an  entirely  different  configuration. 

• The  transition  section  between  the  compressor  chamber  and  the  barrel  should  conform  to  a 
proven  nozzle  design  to  reduce  heat  loads  and  choking 

• The  barrel  length  must  be  sized  to  account  for  momentum  losses  to  the  boundary  layer  (i.e., 
barrel  lengths  beyond  400  diameters  should  be  avoided) 

• The  piston  motion  must  be  tailored  to  achieve  a compromise  between  optimum  projectile 
acceleration  and  tolerable  heat  loads  (e.g.,  a CBPL  operation  for  high  projectile  acceler- 
ation until  midway  in  the  cycle,  and  thereafter  SWL  operation  with  decreasing  base 
pressure  to  limit  heat  loads  on  materials) 

t The  steel  barrel  (and  possibly  the  steel  compressor)  must  be  protected  from  expected  heat 
loads;  in  conjunction  with  a tailored  piston  motion,  liner/inserts  made  of  tungsten  (or 
possibly  an  ablative  material  such  as  carbon)  should  be  evaluated  for  response  to  heat 
loads 
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APPENDIX  A 

EQUATION  OF  STATE  FOR  HYDROGEN 


In  a pressure  range  of  1 kb  to  15  kb  and  a temperature  of  less  than  7,000°K,  hydrogen  is  a 
reasonably  dense  gas  that  is  essentially  diatomic  with  some  dissociation.  Under  these  conditions, 
this  gas  is  characterized  as  a Van  der  Waals  medium  which  is  subject  to  the  effects  of  molecular 
attraction  and -covol urae,  and  influenced  by  vibrational  excitation,  dissociation,  and  ground  state 
electronic  excitation. 


Accordingly,  the  equation  of  state  is: 


pRT  s 2cx+(l-u) 


a(l-g)p 

l-(l-a)pb  RT 


1 


(A-1) 


which  incorporates  the  law  of  partial  pressures  for  hydrogen  dissociation: 


2H 


The  terms  in  Equation  (A-1)  are: 

T = Temperature 

a = degree  of  dissociation 

a = intermolecular  attraction  parameter 

b = molecular  covolume  parameter 

R = gas  constant  per  gram  for  molecular  hydrogen 


In  terms  of  the  liquid-vapor  critical  point  for  molecular  hydrogen: 

V 


27 


; b = ^ 


where  V^  is  the  experimentally  determined  critical  volume 


= 34.57  cmVgm  (Reference  A-1] 


T^  = 33.2°K 


Pj.  = 12.98  X 10  kb 


R = 4.12475  X 10’^  hydrogen 
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The  degree  of  dissociation  of  the  gas  is  derived  by  considering  the  gas  in  equilibrium  with 


Gibbs  free  energy  at  a minimum  (Reference  A-2). 

, exp(-0,/T) 


g 

1-g 


0.0392 


v?{.- 


exp(-0^/T) 


(A-2) 


where  T = Temperature 

= characteristic  dissociational  temperature  = 51  ,965°K 
= characteristic  vibrational  temperature  = 6,324°K 

The  constant  ''0.0392"  is  evaluated  from  fundamental  constants 


\ g^^ 


where  = mass  of  hydrogen  atom  = 1,67  x 10*  gm 
n 

k = Boltzmann  constant  = 1.38  x 10~^^erg/°K 

— 27 

h = Plancks  constant  = 6.625  x 10  erg-sec 

0JJ  = characteristics  rotational  temperature  = 87.5°K  (Reference  A-1) 

U 

g^  = ground  state  electronic  partition  function  for  H-atom  = 2 

U 

gg2  = ground  state  electronic  partition  function  for  H2-molecule  = 1 
The  internal  energy  per  gram  of  hydrogen  is  given  by: 


e = . RT  I 3g.  (l-a)j  5/2  . 


where  Z = 0^/T 


The  Helmholtz  free  energy  per  gram  for  hydrogen  is  given  by: 


t r /zttm  k\' 

f = -2aRT 


4M, 


+ ?,n 


2gp 


^2T-V 


(1-g)  RT  fn 


(A-3) 


(A-^) 


“A.Mnky/^  Mh‘ 

+ tn 

§R_ 

(l-g)o(l-exp(-z)) 

RT  ( : 
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and  is  derived  from  f = X n.f. 

1 1 


where  n^  - mass  fraction 


f.  = -kT 


I Jn;p  +1  I 


IP  = partition  function 

Partition  functions  for  translation,  rotation,  ground  state  electronic  excitation,  vibration, 
dissociation,  covolume,  and  Van  der  Waals  attraction  are  obtained  from  References  {A-1  and  A-2) 

The  entropy  per  gram  of  hydrogen  is  given  by: 


The  sound  speed,  C and  isentropic  index,  are  calculated  numerically  using  (A-1),  (A-2). 


and  (A-5) 


[£]*'■ 

L -IaS  = 0 


AT  = 1°K 


V = £ 

'sent  p 


AS  = 0 
AT  = 1°K 


Figures  (A-1)  through  (A-5)  depict  the  behavior  of  hydrogen  as  defined  by  Equations  (A-1)  through 
(A-7).  Figure  A-1  is  a modified  "Mollier"  diagram  of  specific  internal  energy,  e vs.  scaled  entropy, 

S/R  for  the  range 

2,000°K  < T < 6,000°K 
0.0005  gm/cm^  1 P 1 0-05  gm/cm^ 

The  other  figures  display  pressure,  degree  of  dissociation,  sound  speed  and  isentropic  index  varia- 
tions with  density. 
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APPENDIX  B 


EXAMPLE  OF  CALCULATION  PROCEDURE  FOR  HYBRID  LAUNCHER  BALLISTIC  CYCLES 

Using  the  methods  described  in  Subsections  3.1  through  3.3,  performance  estimates  can  be  made 
for  various  configurations  and  initial  conditions  of  the  PI  launcher.  For  each  estimate,  ideal  driver 
performance,  complete  injection,  and  an  adiabatic  reservoir  are  assumed.  The  launch  tube  length  is 
assumed  to  be  300  bore  diameters,  a length  consistent  with  experimental  and  theoretical  work  on  maxi- 
mum usable  barrel  lengths,  as  determined  by  boundary  layer  losses. 


As  an  example,  a performance  estimate  can  be  made  using  the  initial  conditions  of  PI  Test  7. 

In  this  test  both  the  drivers  and  compressor  section  were  initially  at  350  psia  and  288°K  and  the 
mass  of  hydrogen  in  the  system  was  83  grams.  The  drivers  are  assumed  to  operate  ideally,  generating 
a shocked  gas  pressure  of  0.97  kilobars.  All  of  the  gas  is  assumed  to  be  injected  into  the  compressor 
so  that  the  injection  state  prior  to  compression  is  determined  by  its  density,  ' 0-0067 

gms/cm^  and  internal  energy,  e^  = 0.329  e.u./gm  (1  e.u.  = 10'^  ergs).  From  Equations  (A-1)  to 
(A-5)  the  other  state  variables  prior  to  compression  are: 


P3  = 0.8  kb 
T3  = 2,780°K 
S/R  = 17.25 

The  2,724  gm  (6  lb)  compressor  piston  is  accelerated  to  0.914  km/sec  (3,000  fps)  and  has  a 
total  energy  of  11.4  e.u.  when  it  begins  to  compress  the  gas.  The  total  energy  of  the  gas  upon  in- 
jection is  27  e.u.,  so  the  energy  contribution  of  the  explosive  drivers  is  substantial. 


The  piston  isentropical ly  compresses  the  gas,  as  it  decelerates  to  rest  172  cm  downstream  of 
the  injection  ports  (102  cm  from  the  start  of  the  transition  section).  The  final  compressed  state 
(reservoir  state)  is  characterized  by  its  density  = 0.0159  gm/cm^  and  entropy  S^/R  = S^/R  = 17.25. 
The  other  reservoir  state  variables  are  calculated  from  (A-1)  to  (A-7),  and  are: 


P4  = 2.98  kb 
T4  = 3,715°K 
C4  = 5.34  km/ sec 
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A 


Y4  = 1-518 


= 0.46  e.u./gm 

~ n 

The  muzzle  velocity  is  calculated  for  the  following  conditions: 

Mp  = 86  gms 

L = 300  bore  diameters  = 762  cm 
A = I X 2.54^  = 5.07  cm® 

P4  = 2.98  kb 

^4 

— = 3.52  km/sec 
^4 

G/Mp  = 83/86  ~ 1 

(The  Seigel  curves  for  = 2 
and  Dg.!  = 5 are  nearly  identi- 
cal, so  either  ran  be  used 

= 4.93  km/sec  for  Y4  =1.4 

The  ballistic  parameter,  x = 0.55,  is  used  in  Figure  6 for  G/Mp  = 1 and  gives  a dimensionless 
velocity  U = 0.75.  Thus,  estimated  projectile  velocity  is  iTp  = UC4  = 3.70  km/sec.  This  compares 
with  Pi's  computed  velocity  of  3.65  km/sec  (Reference  I)  determined  by  the  PISCES  code  for  the  same 
initial  conditions. 

The  entire  calculation  can  be  repeated  with  the  covolume  forces  suppressed  to  give  the  fol- 
lowing reservoir  conditions: 

p^  = 2.89  kb 

= 0.0187  gm/cm® 

T^  = 3,71 5°K 
“4  ~ 1% 

C4  = 4.43  km/sec 
-4  = 1-26 
C4/Y4  =3.5  km/sec 

It  is  interesting  to  note  that  P4,  T^  and  C^/y4  are  essentially  unchanged,  while  p^,  and 
Y4  are  considerably  different.  Since  p^  and  C4/Y4  are  essentially  unchanged,  the  calculated  muzzle 
velocity  of  3.67  km/sec  with  covolume  forces  suppressed  is  very  close  to  that  calculated  with  the 
full  hydrogen  equation  of  state,  3.70  km/sec. 
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